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1. - Introduction. 

The relevance of many-body effects, over and above the independent- 
particle approximation, on the optical properties of semiconductors has been 
increasingly recognized over the years. Specifically, the collective effect of 
screening has long been known to be essential to account for the reduction of the 
electron-hole attraction in an exciton and for the polarization accompanying a 
single-particle excitation. Only more recently however, it has been possible to 
give quantitative account of these phenomena by combining Green's function 
techniques of quantum field theory with the description of band structures in 



(*) Present address: Scuola Normale Superiore, 56100 Pisa, Italy. 

1 



2 



G. STRINATI 



terms of localized (molecularlike) orbitals [1-3]. The theory of elementary 
excitations in crystalline semiconductors with electronic states that are to some 
extent localized, which has emerged from these efforts, can be considered rather 
well established by this time. 

Objective of this paper is to give a pedagogical review of the theoretical 
framework underlying the calculations of many-particle properties in semi- 
conductors which are based on the concepts of single- and two-particle ex- 
citations. Biased by the beliefs that, quite generally, only through detailed 
knowledge of the theory one can correctly formulate the problems and that, for 
the specific physical phenomena of interest, the Green's functions method 
provides the very language to describe them, we shall put special emphasis on 
the working procedures that are important for appreciating the subtleties of the 
method, but are not commonly described in the literature. For these reasons, 
results of numerical calculations will be only sketchly presented and references 
will be given to the original papers for additional details. In particular, when 
relating theoretical results with experimental quantities, the direct link between 
the two-particle Green's function and the correlation functions of linear response 
theory will extensively be exploited. It is, in fact, the focus on the correlation 
functions that renders the Green's functions method quite efficient and practical 
by avoiding the calculation of redundant information. 

To keep the presentation as compact as possible, we shall avoid the for- 
mulation of Green's functions theory in terms of conventional diagrammatic 
techniques [4], but we shall rather adopt the alternative formulation in terms of 
functional derivatives techniques that reduces the many-body problem directly 
to the solution of a coupled set of nonlinear integral equations [5]. In this way, 
the screening mechanism which is so important in a semiconductor will be 
introduced from the outset in the formal theory by replacing the ordinary 
Coulomb interaction between electrons by a modified (time-dependent) 
interaction that takes into account the polarization of the medium represented by 
the remaining electrons. 

We confine ourselves in this paper to the zero-temperature limit which is 
appropriate to describe the optical properties of a semiconductor, although 
extention to finite temperatures is feasible [6]. The system we consider are N- 
interacting electrons moving in the static potential of the ions. The spatial 
symmetry of this potential will be left as much as possible unspecified but we will 
restrict eventually to insulting systems with crystalline (space group) 
symmetry, while actual calculations will be presented for cubic covalent 
semiconductors. No coupling with phonons as well as no magnetic and relativistic 
effects will be considered throughout. 

This paper is organized in two parts. The first part (sect. 2-7 and appendices 
A-E) introduces the theoretical tools and the associated calculation methods and 
approximation procedures which connect the study of the optical properties of 
semiconductors to the Green's functions method. The second part (sect. 8-11 and 
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appendices F-H) deals more specifically with the computation of quantities (such 
as the optical absorption coefficient) characterizing single- and two-particle 
excitations in (covalent) semiconductors. 

2. - Generalized Green's functions. 

We consider a (nonrelativistic) iV-eleetron system which interacts with a 
static potential V(r). The corresponding second-quantized Hamiltonian has the 
form 

(2. 1) H = j dx f f (x) h(r) f{x) + || dx dx' f\x) f\x') v(r, r') f(x') Y(x) . 

In eq. (2.1) x signifies the set of space (r) and spin (£) variables, the Y(x) are field 
operators, h(r) is the one-electron Hamiltonian 

(2.2) Mr) = _^L V 2 + F(r) 

2m 

and 

(2.3) v(r>r > )= *^ 

\r— r | 

is the Coulomb interaction [7]. 

We consider also an external (scalar) potential U(x, x'; t) which is local in 
time but nonlocal in space and which couples bi-linearly to the field operators. 
The interaction Hamiltonian is thus taken of the form 

(2.4) H'(t) = j dx dx' Y\x) U(x, x'; t) ¥(x') , 

which is guaranteed to be Hermitian by requiring U(x, x'; t) to be Hermitian in 
the x variables at any t. This potential is further assumed to vanish as 00 . 

The introduction of the interaction (2.4) may be considered as a purely formal 
tool to generate in a compact form the equations of motion for the Green's 
functions together with a number of useful relations. It will be understood, in 
fact, that the external potential will be allowed to vanish at the end of the 
calculation. Nevertheless, the formalism we shall develop could be used as well 
to follow the time development of the Green's functions of the system under the 
action of a physical external (local) potential of the form 

(2.5) U(x, x'; t) = U(x, t) S(x, x') . 

More general forms of coupling to scalar and vector (electromagnetic) potentials 
will be considered in sect. 5 and in appendix C. 

The interaction Hamiltonian (2.4) enables us to introduce an interaction 
picture by defining the time dependence of the field operators with respect to the 
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unperturbed Hamiltonian (2.1): 

(2.6) f(l) = ?"(*i , <i) = exp [iHt! Ih] Y(x{) exp [ - iHtj. lh\ . 
Similarly, we write 

(2.7) H[{t) = exp [iHt/h] H'(t) exp [- iHt/h] = 

= jdxdjc' f f (x, t + ) U(x, x'; t) f(x', t) , 

where for later convenience t + stands for t + 8 + ), and we introduce the 
(formal) operator 

(2.8) S = exp|-| } df#KoJ. 

The generalized single- and two-particle Green's functions are then defined 
to be 

(2.9) GA2)= _l(NmSm)t'(MN) 



h <JV|2TS]|JV> 



and 



(2.10, «!, 2; r , *>-| - ft <mSnl) r<2) T ' m rm]lN) 



h) (N\T[S\\N) 

respectively, where \N) denotes the ground state of the unperturbed iV-electron 
system and T is Wick's time-ordering operator which includes a minus sign for 
any permutation of (fermion) field operators. Equations (2.9) and (2.10) reduce to 
the definitions of the ordinary single- and two-particle Green's functions [4] 
when the external potential U is allowed to vanish. Notice also that all the U- 
dependence in the generalized Green's functions is contained in the operators S. 



3. - Equations of motion for the generalized single-particle Green's function. 

The presence of the time-ordering operator in eqs. (2.9) and (2.10) requires 
us to consider, besides the S-matrix T[S], the time evolution operator in the 
interaction picture T[S(t a , t b )], where 

(3. 1) S(t a , t b ) = exp j - 1 | dtHKt) 1 



(with the understanding that the operator (3. 1) makes sense only within a time- 
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ordered product). This operator satisfies the following relevant properties: 

(3.2) i) T[S(t a , Q] = T[S(t a , t b )]T[S(t b , y] 
(group property); 

(3.3) ii) T[S(t a , t b ) F(l) F(2)] = T[S(t a , «] F(l) F(2) T\$(k, 4)] , 
for t a >t 1 >t 2 >t b ; 



(3.4a) 



(3.46) 



iii) 



4~ T lS(t a , t b y\ = - lmt a ) T[S(t a , t b )] , 
at h 



A- T[S(t a> = i T[§(t a , t b )] H[{t b ) . 
at b n 



The equation of motion for Gi(l, 2) and its adjoint are then obtained by 
taking the derivative of Gi (1, 2), alternatively, with respect to ti and t?,. After 
the time dependence of the various factors entering the definition (2.9) is made 
explicit, in addition to eqs. (3.1)-(3.4) we require the equations of motion of the 
field operators 



(3.5a) 



^-F(l) = -f 
all h 



Ml) + |d3v(l, 3)F + (3)F(3) 



(3.56) -Ltf* ( 2) = i 

3t2 n 



h{2) F + (2) + F f (2) \ d3v(2, 3) F*(3) F(3) 



as well as the identity 
(3.6) 



where 6 is the unit step function. In eqs. (3.5) we have introduced the notation 
v(l, 2) = v(xi, r 2 )S(t 1 — ti) and h(l) = fefo), while the integrals extend over space-, 
spin, and time variables which are collectively denoted by 1, .... 
Manipulations lead then to the following equations: 



(3.7a) 



(3.76) 



ih-^-hO.) dd, 2)-/ d3tf(l, 3)G X (3, 2) 



+ 



+ ihj d3v(l, 3) G 2 (l, 3 + ; 2, 3 ++ ) = 2) , 



- tft^- - m Gi(l, 2) - / d3G 1 (l, 3) [7(3, 2) 



+ 



+ ift|d3-y(2, 3)G 2 (1, 3~; 2, 3") = <5(1, 2), 
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where 3* implies that the time variable t 3 is augmented (diminished) by a 
positive infinitesimal, and we have introduced the notation 



(3.8) 



17(1, 2) = U(x u x 2 ; tOM-fe). 



A characteristic feature of eqs. (3.7) is that, besides the single-particle 
Green's function, they involve also the two-particle Green's function. A whole 
hierarchy of equations involving higher-order Green's functions can thus be 
generated in this way [8]. It is customary to replace this hierarchy of equations 
by a coupled nonlinear set of integro-differential equations connecting the single- 
particle Green's function to the self-energy operator (£), the (irreducible scalar) 
vertex function (T), and other derived quantities like the (irreducible) polar- 
izability (%), the dynamically screened interaction (W), etc. To this end, we begin 
by eliminating formally the two-particle Green's function from eqs. (3.7) by 
utilizing the following functional derivative identity: 



(3.9) 



G 2 (l, 3; 2, 3 + ) = G 1 (l, 2)^(3, 3 + ) ■ 



SGid, 2) 
8*7(3) 



which holds in this form for generalized Green's functions defined in the presence 
of the local potential U(l) (cf. eq. (2.5)). Potentials of the more general form (3.8) 
will instead be needed to generate the equation of motion for the two-particle 
Green's function. Equation (3.9) is derived in appendix A where the calculus with 
functional derivatives is also briefly reviewed. 

Inserting the identity (3.9) into the equations of motion (3.7) yields two 
functional differential equations: 



(3.10a) 

and 
(3.106) 



Gid, 2) - 
SGjd, 2) 



-ih\a3v(v, 3) J =$a, 2) 



~ ih ~ik~ h{2) ~ U{2) +ih i d3v(2 ' 3) Gi(3 ~' 3) 

-ih fd3u(2-, 3) 



8*7(3) 

G.a, 2)- 

8Gi(l, 2) 
S£/(3) 



= *(1, 2). 



Equations (3. 10) are not yet in a form suitable for taking the limit as ?/— > 0, since 
this process would require explicit knowledge of the functional dependence of Gi 
on U itself. To avoid any reference to the external potential, it is convenient to 
rewrite eqs. (3. 10) in terms of the self-energy operator which is defined by the 
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two alternative expressions 



(3.11a) 2(1, 2) = 2 H Q, 2) + ih Jd34v(l + , 3) 



8Gi(l, 4) 



(3.11b) 2(1, 2) = 2 H Q, 2) + ift Jd34Gr 1 (l, 4) 



817(3) 
«Gi(4, 2) 



Gr 1 ^, 2), 
v(3, 2"). 



817(3) 

Here 2 H stands for the Hartree contribution to the self-energy 



(3.12) 



2 H Q, 2) = 8(1, 2) 



■ift Jd3v(l, 3)Gi(3, 3 + ) 



while the rest is sometimes referred to as the mass operator and indicated by 
M(M). The inverse Gf 1 of the generalized single-particle Green's function (2.9) is 
defined by 

(3.13) j dSGtHX, 3)G,(3, 2) = J dSG^l, 3)Gr x (3, 2) = 8(1, 2). 

Demanding the left and right inverse of Gi to be equal requires the two 
expressions (3.11) to coincide. This property can be verified for the exact self- 
energy through a diagrammatic expansion; its relevance will appear when it will 
be imposed as a constraint on the choice of the approximate form of the self- 
energy (sect. 6). 

Equations (3.10) can now be cast in the form 



(3.14a) 



ih-Z—KD-UQ) 
ati 



or, alternatively, 

9 



(3.146) 



ih^--h(2)-U(2) 



Gi(l, 2)- Jd3£(l, 3)GA 2) = 8(1, 2) 



Gid, 2) - 1 d3Gi(l, 3)2 (3, 2) = 8(1, 2) . 



Equations (3.14) express Dyson's equation in differential form (cf. appendix B) 
and give to X the physical meaning of a nonlocal and energy-dependent effective 
single-particle potential, as it will be discussed in sect. 9. Hereafter we shall limit 
to consider only one of the two equations (3.14), since the other is implied by the 
condition 2 = 2. 

We proceed now to eliminate any explicit reference to the external potential 
in the self-energy operator (3.11). To this end, we introduce the following 
auxiliary quantities. 

i) The total (classical) potential V is defined as 
(3. 15) V(l) = U(l) -ihj d3v(l, 3) Gi(3, 3 + ) . 
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Regarding alternatively Gi as a functional of V instead of U and using the 
«chain rule» (A.8) as well as eq. (A. 7) yield 

(3.16) M(l, 2) = ~ihj dS45v(l + , 3)^(1, 4) ^ j^. 

ii) Equation (3.16) suggests to introduce a scalar (irreducible) vertex 
function f defined as 

( 3.1 7 ) P (1)2; 3),-i™ = ,(M),(2,3) + SM(1 ' 2) 



where the right side has been obtained by combining eqs. (B.4), (3.11), and 
(3.15). Since the mass operator M depends on the external potential U (or 
alternatively on V) only through its dependence on the generalized single- 
particle Green's function (cf. eqs. (3.11)), we can again use the «chain rule» for 
functional differentiation together with eq. (A. 7) to obtain an integral equation 
for f: 

(3.18) ni,2;3) = ,(l,3),(2,3) + /d45^« = 

= &{!, 3)*(2, 3)+ f d4567 ^. ( . 1 'f. ) G 1 (4 > 6)^(7, 5)r(6, 7; 3). 

J SCri(4, 5) 

This can be considered as the equation defining f, whereby the limit U—> can 
explicitly be taken. 

iii) Equation (3.16) suggests to introduce also an inverse (longitudinal) 
dielectric matrix e" 1 defined as 

(3.19) ^^-m>- 

From the definition (3.15) and from the identity 

taw -« wl ,n.fff. W) ), 

(N\T[S]\N) 

where p(l) = F + (l) T(l) is the density operator, we can express 

(3.21) 6" HI, 2) = *(1, 2) + jd3v(l, 3)^|p 

The functional derivative of the average density with respect to the external 
potential in eq. (3.21) is a measure of the polarizability of the system. Its 
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connection with the corresponding quantity of linear response theory (cf. sect. 5) 
is most readily drawn by utilizing the identity (3.9): 



(3.22) 



X (l, 2) = ^^-=ifc[G 2 (l, 2; 1 + , 2+)-G 1 (l, l + )G 1 &, 2 + )] = 



SU(2) 



i (N\T[Sp'a)p'(2)]\N) 
h 



{N\T\S\\N) 



where we have introduced the density deviation operator 
(3.23) p'(l) = p(l)-(p(l)>. 

From the last line of eq. (3.22) x is seen to be symmetric under the 
interchange of its arguments. Moreover, it is convenient to single out from x the 
part which is irreducible with respect to the bare Coulomb potential v by 
regarding the average density as a functional of the total potential V: 



(3.24) x (l, 2) = /d3 



= /d3 



8<p(l)> svm 
87(3) 8t7(2) 



-J 



d3 



8<P(D) 



87(3) 



e-H3, 2) = 



8<p(4)>' 



*(3, 2) + \ d4v(3, 4)- s[7(2) 

= z(l, 2) + / d34x(l, 3)tK3, 4) z (4, 2) , 
where we have defined the irreducible polarizability x to be 

«<#!)> 



(3.25) 



X(l, 2)i 



SV(2) 



Equation (3.24) can be regarded as an integral equation to be solved for x once 
the kernel x is specified. Knowledge of x, in turn, can be related to that of the 
vertex function (3.17) through the identity (A. 7) 



(3.26) x0-,2) = -ih 



SGid, 1 + ) 
87(2) 



= ift Jd34Gi(l, 3) 



SGlKS, 4) 



Gi(4, l + ) = 



87(2) 

= -ihj d34Gi(l, 3)Gi(4, l)r(3, 4; 2). 



Notice that x also is symmetric under interchange of its arguments, and that we 
can express the {longitudinal) dielectric matrix e in terms of x as follows: 



(3.27) 6(1,2) = 



817(1) _ 8 
87(2) 87(2) 



7(l)-/d3v(l, 3) (p(3))J = 

= <?(!, 2)-Jd3t*l, 3)z(3, 2), 
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where 

(3.28) | d3e~ \1, 3) e(3, 2) = J d3e(l, 3) e~K3, 2) = <y(l, 2) . 

Equation (3.28) follows from eq. (A.3) and the «chain rule» (A.8). 

iv) It is further convenient to introduce the dynamically screened 
interaction W defined as 

(3.29) W(l, 2) = |d3e- 1 (l, 3)v(3, 2). 
We can express alternatively 

(3.30) W(l, 2) = v(l, 2) + | d34«(l, 3) #3, 4) W(4, 2) = 

= v(l, 2) + Jd34v(l, 3) z (3, 4)v(4, 2), 

where the first line can be interpreted as an integral equation defining W, while 
the last line represents its formal solution. From the last line it also follows that 
W(2, 1) = W(l, 2). 

With the definitions (3.17), (3.19), and (3.29) the mass operator (3.16) can be 
cast in its final form 

(3.31) M(l, 2) = ih \ d34W(l + , 3) Gjd, 4)2X4, 2; 3) . 

Equations (3.14), (3.18), (3.26), (3.30), and (3.31) constitute the set of coupled 
equations which link M, f, x, and W. They are represented graphically in fig. 1 
together with the expression for M (cf. equation (3.11&)) 

(3.32) M(l, 2) = ihj d34f(l, 4; 3) Gi(4, 2) W&, 2") . 

Notice that the limit f/-»0 can be readily taken in these equations. 

All quantities considered thus far are still exact. Approximations can be 
generated either by expressing the set of coupled equations as (infinite) series in 
terms of Gi 0) and v, thereby reproducing the Feynman-Dyson perturbation 
theory [4], or by truncating the set of coupled equations by making specific 
ansatz on the functional form of the mass operator M in terms of the self- 
consistent G 1 and W [5, 6, 8]. Different approximations will be relevant to 
different physical situations, and for any given approximation it will be 
important to check whether it satisfies rather general conservation criteria such 
as the conservation of the number of particles, the total energy, etc. A 
prescription to generate such conserving approximations will be discussed in 
sect. 6. 
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G, (0) (1,2)= 



M (1,2) = 



M(l,2) = 



I'd, 2) = • • 



1 2 
W(l,2) = rwvww^ 



r(i,2;3)= W/w>» 



5MQ,2) 
66,(3,4) : 



*0,2) = 




a) 



2 1 



2 1 
-• + •- 



=» + 




1 2 1 2 1 

•x/vvvvxrty"* = • • -f- • 




■v2 1 




k2 1 




Fig. 1. - a) Graphical symbols corresponding to the relevant quantities needed to 
represent b) the coupled set of equations for M, f, x, and W. 



4. - Bethe-Salpeter equation for the two-particle Green's function. 

The equation of motion for the two-particle Green's function, which is the 
analog of the Dyson's equation (B.6) for the single-particle Green's function, can 
be most readily obtained by introducing the two-particle correlation function 
defined as 



(4.1) L(l, x' t; 2, xt + ) = - G 2 (l, jc' t; 2, xt + ) + Gi(l, 2) Gib' t, xt + ) . 
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In fact, for generalized Green's functions evolving in the presence of an external 
potential of the type (3.8), L can be expressed by means of the identity (A. 12), 
namely, 

SGid, 2) 

(4.2) L(l, x' t; 2, xt + ) = 



SU(x, x';t)' 

A combination of eq. (4.2) with eqs. (A. 7), (A.8), and (B.4) then yields 

r SC~ 1 (3 4) 

(4.3) L(l, x' t; 2, xt) = - J d34G!(l, 3) ^ ' , ' G,(4, 2) = 



SU(x, x'; t) 

$2(3, 4) 

*<& - U) ${U ~ t) S(x s , x) S(x 4 , x') I §U{x x ,. t) 



G l (4, 2) = 



= Jd34G 1 (l, 3) 

= Gi(l, xt) Gj(x' t,2) + j d3456Gj(l, 3) Gi(4, 2) 3 ^ ^ L(6, 5, xi + ) . 

This is an integral equation for L, whereby x, x', and t play the role of external 
variables and the kernel 

5Z(3, 4) 

(4.4) 5(3, 5; 4, 6)=- 



SGtf, 5) 



represents an effective two-particle interaction. Notice that the limit as U—> 
can be explicitly taken in eq. (4.3). Notice also the following properties. 

i) The topology of its diagrammatic structure (cf. fig. 2a)) implies that eq. 
(4.3) can be generalized to hold for arbitrary values of the external time variables 
t and V and not just in the limit t' = t~. We can then write compactly 

(4.5) L(l, 2; 1', 2') = G 1 (l, 20^(2, 1') + 

+ J d3456G!(l, 3)Gx(4, l')S(3, 5; 4, 6)L(6, 2; 5, 2'). 

Equation (4.5) is known as the Bethe-Selpeter equation for L [8]. 

ii) It is convenient to single out from the outset the Coulomb term in the 
effective interaction E by breaking up the self-energy 2 as the sum of the 
Hartree term (3.12) and the mass operator M (cf. eqs. (3.11)), to obtain 

s 5M(3, 4) 

(4.6) J?(3, 5; 4, 6) = - ihS(3, 4)<?(5, 6)v(3, 6) + ^ g) ■ 

This splitting will enable us to identify the irreducible (or proper) part of L as 
well as other correlation functions (appendix C). 
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6G, 



2' 1 , l ' 1 
= + 



2 l' 2 1' 



2' 1 



2 1' 4 
1 3 

+ 



1' 



4 

a) 



6) 



6M 
66, 





2' 1 



2 1' 



2' 1 



2 1' 

c) 



Fig. 2. - Bethe-Salpeter equation for a) the two-particle correlation function L, b) the 
bound part of the two-particle Green's function SG 2 , and c) the many-particle ^-matrix. 



iii) A connection between the integral equations (4.5) for L and (3.18) for r 
can be drawn from the identity 

(4.7) L(X, 3; 2, 3 + ) = / d45Gx(l, 4) G 1 (5, 2)r(4, 5; 3) , 

where we have introduced the scalar (reducible) vertex function r defined as 

(4.8) 7X1, 2; 3) s J d4T(l, 2; 4) e~\A, 3) = 

s gtHi, 2) gy(4) _ sGiKi, 2) 
«v(4) ms) sum ' 



= - Jd4- 
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The last line of eq. (4.8) has been obtained from the definitions (3.17) and (3.19) 
and the «chain rule» (A.8). 

Equation (4.5) is not the only form of the Bethe-Salpeter equation one finds 
in the literature. Two alternative forms are shown graphically in fig. 2b) and 2c). 
They can be obtained from eq. (4.5) as follows. The two-particle Green's function 
G 2 can be separated into a free part corresponding to the propagation of two 
single-particle excitations totally independent of one another and a bound part 
SG 2 [9] 

(4.9) G 2 (l, 2; 1', 2') = ^(1, l')C(2, 2') - G,(l, 2')G 1 (2, 1') + SG 2 (1, 2; 1', 2') , 
whereby the bound part §G 2 satisfies the integral equation [9, 10] 

(4.10) 8G 2 (1, 2; 1', 2') = - \ d3456G 1 (l, 3)^(4, V)E{Z, 5; 4, 6) • 

•^(6, 2')G!(2, 5) + | d3456G!(l, 3)^(4, V)E(S, 5; 4, 6)SG 2 (6, 2; 5, 2'). 

Otherwise, one can introduce the many -particle T-matrix defined as the solution 
of the integral equation [9, 11] 

(4.11) 7X1, 2; 1', 2') = S(1, 2; 1', 2') + 

+ Jd3456£(l, 4; 1', 3)Gj(3, 6)Gi(5, 4)7X6, 2; 5, 2'). 

Solving eq. (4.11) is equivalent to solving eq. (4.10) since the bound part SG 2 can 
be obtained by attaching a single-particle Green's function to each end point of T: 

(4.12) SG 2 (1, 2; 1', 2') = 

= - / d3456G 1 (l, 3)^(4, l')7X3, 5; 4, 6)^(6, 2')G 1 (2, 5). 

The spin variables can be eliminated from explicit consideration whenever 
the Hamiltonian is spin independent. The procedure is trivial for the single- 
particle Green's function, the vertex function, and the polarizability, but it 
requires some care for the two-particle Green's function or the 7 1 -matrix. In 
particular, the effective interaction entering the Bethe-Selpeter equation for the 
singlet channel differs from that for the triplet channels (appendix D). 

Comparison of eqs. (3.22) and (4.1) shows that the polarizability x can be 
obtained as a degenerate form of the two-particle correlation function L, namely, 

(4.13) x (l, 2)=-tftL(l, 2; 1 + , 2 + ). 

In this limit the Bethe-Selpeter equation (4.5) is equivalent to the set of 
equations (3.24), (3.26), and (3.18), thereby providing information on the density 
fluctuations of the system and related physical quantities (such as plasmons). 
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The full equation (4.5), on the other hand, is relevant to the study of bound-state 
problems (such as excitons) in the spirit of the original work of Gell-Mann and 
Low [12] (cf. sect. 11). 

5. - Connection with linear-response theory. 

The time-ordered correlation functions introduced thus far (cf. appendix C) 
do not directly represent the effects of the coupling between a system of 
particles and an external agent when causal boundary conditions are required. 
Specifically, if the system was in the ground state before the external agent 
acted, the linear response of the ground-state expectation value of a given 
operator 6{x, t) is known to be given by a Kubo formula [13] 

t 

(5.1) $(6(x, t))=± \ dt'{N\[H{(t'), 6i(x, t)]\N) . 

Here 

(5.2a) H[{t) = exp [iHt/h] H'(t) exp [- iHt/h] , 

(5.26) 6i(jc, t) = exp [iHt/h] 6(x, t) exp [- iHt/h] , 

where H'(t) is the (weak) time-dependent interaction Hamiltonian between the 
system and the external agent. In particular, for an electromagnetic field H'(t) is 
given by eq. (C.l). In this case one is interested in the linear response of the 
density and current operators, given by eq. (C.2) and by 

(5.3) /(x, t) =j(x) - —A{r, t)p(x) , 

mc 

respectively (cf. eq. (C.3)). To first order in the fields, one obtains for the total 
density and current: 

(5.4) (p(l)) A , ? = (N\p(l)\N) + q / d2^ R (l, 2) 9 (2) - 1 x E (l, 2) • A(2)j , 

(5.5) (JO)) a,, = ~~ (N\p(X)\N) A(l) + 

+ qj d2^ E (l, 2) 9 (2)-±x*(l, 2)-A(2)j, 

where we have introduced the retarded counterparts of the correlation functions 



16 



G. STRINATI 



(C.23) and (C.24): 



(5.6a) ZR (l, 2) = - j- (N\[p\l), p'(2)]\N) 6(t x - t 2 ) , 



(5.66) %(1, 2) = - j- (N\[j'(l), p'(2)]\N) 6(t x - t 2 ) , 



(5.6c) 2 R (i, 2) = - j- (N\[j'(l), f(2)]\N) 6(t x - t 2 ) . 

Here the unit step function 6(t x - t 2 ) enforces the casual boundary conditions 
as it limits the knowledge of the interaction between the system and the external 
agent to antecedent times; the time-ordered correlation functions (C.23) and 
(C.24), on the other hand, require also the knowledge of the future of this 
interaction and thus do not have bearing on the experimental situation. 
However, the retarded correlation functions cannot be calculated through the 
set of coupled integral equations developed in appendix C (or, alternatively, by 
the Feynman-Dyson perturbation series) because the identities (C.4) and (C.5) 
hold only for time-ordered products of operators (or, alternatively, Wick's 
theorem applies only to these products). A connection between the two sets of 
correlation functions is then in order. This connection can be drawn by looking at 
the respective Lehmann representations after Fourier tranforming the time 
dependence. 

For a general pair of time-ordered and retarded correlation functions 

(5.7) C^(l, 2) = -±-(N\T[A'(l)B'(2)]\N) , 

n 



(5.8) Ci B (l, 2) = -^ (N\[A'(l), B'(2)]\N) 6(t, - t 2 ) , 

where A and B are Hermitian Bose-like operators, we obtain 

A s ( Xl )Bf(x 2 ) A*(x 1 )B,(x2) 



(5.9) C^(x u x 2 ; co)=Z 



\ ha, + E - E s + iS hoj + E s - E - iS 



tab,„ ... ^ ^/ A s (x x )B*(x 2 ) Af(x 1 )B s (x 2 ) 



(5.10) C^Xu x 2 ; w ) = 2 



* \fiu + E -E S + iS hoj + E S ~E + iS 

Here the index s labels the eigenstates \N, s) of the unperturbed Hamiltonian H 
(s = corresponds to the ground state \N)), (E s -E )>0 is an excitation energy, 
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S-+0 + , and 



(5.11) 



A s (x)=(N\A(x)\N, s). 



Note the following properties (for real w): 



(5.12) 



i) Cr £ (x 1 , x 2 , —oj) — C^iXi, x 2 , co) 



which implies that its real part is an even and its imaginary part is an odd 
function of w; 

(5. 13) ii) CT(x b x 2 ; - w ) = C$ A {x 2 , Xl ; w) ; 

(5.14) iii) Ci B (x 1 ,x 2 ; ( o) = Cf s (x 1 ,x 2 ;<o), 

for &j > 0. One can thus evaluate Cf^Xi, jc 2 ; a>) for w > by standard many-body 
techniques and then obtain Cr^jci, x 2 ; &>) for all values of w through eqs. (5.12) 
and (5.14). 

Further properties can be obtained when the operators A and B coincide 
with either p or j. In this case ^t(*i> x 2 ; w) and %t(*i> *2,' &>) are even functions of 
w, while xt(*i, w) is an odd function of w (provided no magnetic field is 
present). 

6. - Conserving approximations. 

An approximate solution to the coupled set of integro-differential equations 
discussed in sects. 3 and 4 may or may not be consistent with the general 
(number, momentum, and energy) conservation laws satisfied by the exact 
solution. A particular approximation is then said to be conserving if it satisfies 
the restrictions imposed by the conservation laws. These restrictions have been 
first formulated by Baym and Kadanoff [11] in terms of sufficient conditions to 
be satisfied by the approximate self-energy operator U, and successively recast 
in a simple diagrammatic form by Baym [14]. Fulfillment of conservation criteria 
turns out to be essential to correctly describe transport phenomena (as, for 
instance, in the context of the Landau theory of Fermi liquids [9]); however, it 
may be less important in other contexts such as the determination of single- 
particle and bound electron-hole pairs excited states of a crystal. In any case, one 
should be aware of possible violations of conservation laws although it might be 
difficult to fulfill them in practice. 

A sufficient condition to ensure fulfillment of conservation laws is that any 
approximation for the self-energy operator E be «<?-derivable» in the sense that 
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there exists a functional $ such that [14] 



(6.1) 



£<1,2)=- 



o<Z> 



3<?i(2, 1) ' 



$ is meant to be a functional of the bare Coulomb potential v(l, 2) and of the 
generalized single-particle Green's function G^l, 2) (in the presence of an 
external field U), which has thus to be determined self-consistently via Dyson's 
equation (3.14). Examples that show how $ can be constructed to reproduce 
known approximate expressions of 27 are reported in fig. 3. Notice there that 

* 2" = (8*/5G.) 





a) 




b) 





d) 







e) 




+ o 



Fig. 3. - Examples of correspondence between diagrams for <P and 27. Straight lines here 
denote the self-consistent single-particle Green's function Gi and broken lines denote the 
bare Coulomb potential v. 



condition (6.1) requires certain diagrams for E to be taken together, as for the 
terms c) and /). This feature, in turn, implies that the two alternative 
expressions (3.31) and (3.32) for the mass operator M coincide for the given 
approximation, as anticipated in the discussion of eq. (3.11). 
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Conservation laws for the ground state follow directly when condition (6. 1) is 
satisfied. In addition, conservation laws are also fulfilled within linear response 
if, for the given choice of 2, the two-particle correlation function L satisfies the 
Bethe-Salpeter equation (4.5) with kernel S = &2/8Gi (that is, with no terms 
omitted) [14]. In this case we can express (cf. eqs. (4.4) and (6.1)): 



(6.2) 



S(l, 2; 1', 2')=- 



8 2 $ 



= 3(2, 1;2', 1'), 



5Gi(2', 2)8G 1 (1', 1) 

which leads to the symmetry condition 

(6.3) L(l, 2; 1', 2') = L(2, 1; 2', 1') 

for the approximate two-particle correlation function. 

As an example, we consider the number conservation law expressed in terms 
of G 



(6.4) —ihG^l, l + ) + V r 



h 



2im 



(V 1 -V 2 )i/iG 1 (l, 2) 



2=1+ 



A + Aj iftG . (1>2) + ^. (v? _ vi)Gl(1>2 ) 



2=1+ 



To recover this equation for the approximate G x , we begin by subtracting the 
two alternative forms of Dyson's equation (3.14a) and (3.146), to obtain 

(6 - 5) (l: + ^)^ (1 ' 2)+ £ (V - vi)G ^ 2) - 

- (Vir,) + £7(1) - V(r 2 ) - £7(2)) dd, 2) + 

+ / d3(G 1 (l, 3) £(3, 2) -27(1, 3)GA 2)) = 0, 

and then take the limit as 2— » 1 + . In this limit, one can readily show from eqs. 
(3.11) that 



(6.6) 



Urn f dS(G 1 (l, 3)2X3, 2) -27(1, S)G 1 (S, 2)) = 0, 

2^1 + J 



since 2 = 2 for a conserving approximation. Equation (6.4) thus follows. 

Equation (6.6) can alternatively be obtained by noticing that the functional <P 
is, by construction, invariant under the transformation 

(6.7) Gxd, 2)-> exp [iA(l)] G 1 (l, 2) exp [- iA(2)] - 

-Gid, 2) + i(A(l)-A(2))G 1 (l, 2), 
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to lowest order in the small function A. Variation of $ then gives 



(6.8) 



S<2> = 




SGid, 2) 



SG.il, 2 + ) = 



= Jdl22;(2, l)i(A(l)-A(2 + ))G 1 (l, 2 + ) = 

= ij dLl(l + ) /d2(G!(l, 2)2(2, l + )-2{l, 2)d(2, 1 + )) = 



where account has been taken of eqs. (6.1) and (6.7). (The need to replace 2—* 2 + 
in the first line of eq. (6.8) originates from the Hartree-Fock term.) The 
arbitrariness of A yields eventually eq. (6.6). 

Physically, eq. (6.4) reduces to the continuity equation for the ground state 
in the limit of zero external field U [15]. However, eq. (6.4) holds more generally 
for arbitrary U and thus it contains information about the induced density and 
current within linear response. Specifically, taking the functional derivative of 
eq. (6.4) with respect to C7(3) and recalling eqs. (A. 13) and (C.26), yields 



where the limit for vanishing U is now understood. 

Equation (6.9) holds for the approximate correlation functions, provided L 
satisfies the Bethe-Salpeter equation (4.5) with kernel taken according to eq. 
(6.2). To interpret it as the continuity equation for the induced density and 
current within linear response, the connection established in sect. 5 between 
time-ordered and retarded correlation functions has to be recalled. In particular, 
by Fourier transforming the time dependence of eq. (6.9) and making use of the 
properties (5.12) and (5.14), it follows that (cf. eqs. (5.4) and (5.5)) 



(6.10) ~(pa)) u + V 1 -{ja))u= WJ-XRd, 2) + V 1 -z R (l,2)) r U(2) = 0. 



Equation (6.9) can also be obtained as a particular limit of a more general 
equation relating 2 to L which is satisfied by a «$-derivable» approximation, as 
discussed in the next section. 



7. - Gauge invariance, current conservation, and the Ward identities. 



(6.9) -i-ihLil, 3; 1 + , 2, + ) + V l -ih\j!-^ 1 -V 2 )La, 3; 2, 3 + ) 






In this section we derive an identity relating the exact single- and two- 
particle Green's functions as well as any conserving approximation to these 
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functions. This identity (which will be referred to as the «generalized Ward 
identity*) will be shown to follow from the gauge invariance of the theory and to 
reduce to the continuity equation (6.9) for the induced density and current in a 
particular limit. The ordinary Ward identities [9] will also be recovered as 
special forms of the generalized Ward identity. 

To derive the relationship between single- and two-particle Green's 
functions, we begin by coupling the system to the electromagnetic field 
according to eq. (C.l), where we replace [14] 



(7.1) 



A(r, t)->VA(r, t) , 



9 (r, *)^-I|-A(r, t). 
c at 



The resulting interaction Hamiltonian, which corresponds to vanishing electric 
and magnetic fields (gauge Hamiltonian), can then be utilized to define 
generalized Green's functions according to eqs. (2.9) and (2.10). In particular, the 
single-particle Green's function satisfies the following equation of motion: 



(7.2) 



•» 3 , ft 



ih ir + £t (Vi + Vii(D) 2 - V(rd + ih^-A{l) 



3ii 2m 



8<! 



Gxd, 2; A) — 



■j&za, 3;i)Gi(3, 2; A) = S(l, 2), 



where 
(7.3) 



A(l) = ^-A(l). 



Equation (7.2) can be obtained from eq. (3.14a) with the usual replacement in the 
presence of an electromagnetic field 



(7.4) 



V^V + ^-A. 
the 



where now A and 9 are taken according to eq. (7.1). The solution of eq. (7.2) can 
be expressed as [14] 



(7.5) 



Gj.il, 2; A) = exp [- A(l)] Gi(l, 2; A = 0) exp [A(2)] , 



which can be verified at every order of a diagrammatic expansion of 2 in terms 
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of Gi. Equation (7.5) implies that 



(7.6) 



SGid, 2; A) 



3,1(3) 



= (6X2, 3) -<?(!, 3)^(1, 2) 



.1 = 



for the exact generalized single-particle Green's function as well for any 
approximation to it. 

The desired connection with the two-particle Green's function can then be 
established by expressing the functional derivative in eq. (7.6) in an alternative 
form by recalling eqs. (C.4) and (C.5). We write 



(7.7) SGiQ, 2; A) = -ih |d3 3; 2, 3 + )^A(3) + 



+ 



-^-(Vj-VjOLd, 3; 2, 3') 

2im 



3' = 3 + 



•V 3 A(3) 



for an infinitesimal A. Integration by parts in eq. (7.7) yields 



(7.8) 



5GiQ, 2; A) 

oA(3) 



= ih \~rUl, 3; 2, 3 + ) + 

6t s 



+ V, 



2im 



(V 3 -V 3 .)L(1, 3; 2, 3') 



and comparison with eq. (7.6) leads eventually the generalized Ward identity in 
the form 



(7.9) -2-L(l, 3; 2, 3 + ) + V 3 ' 

9< 3 



_a-(V s -V s .)L(l,3; 2, 3') 

2im 



3' =3 + 



■ Uiil, 3)-<y(2, ^Gid, 2), 



where the limit A = is understood. 

This equation holds for any approximate two-particle correlation function L, 
provided it satisfies the Bethe-Salpeter equation (4.5) with kernel E obtained 
from the approximate single-particle Green's function Gi according to the 
prescription (4.4). However, in order to recover from eq. (7.9) the continuity 
equation (6.9) for the induced density and current within linear response, the 
kernel E has in addition to satisfy the symmetry requirement (6.2) which is met 
by any «$-derivable» approximation. This is because the resulting symmetry 
property (6.3) for L is required to relate the restrictions of local charge 
conservation (6.9) and of gauge invariance (7.9) [16]. 
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The identity (7.9) can be obtained in an alternative way by expressing is left- 
hand side in terms of the density and current deviation operators (C.25) 



(7.10) -£-L(l, 3; 2, 3 + ) + V 3 - 



h 



2im 



(Vs-V 8 .)L(1, 3; 2, 3') 



3' = 3+ 



N 



(3) + V 3 -j'(3)W(l) 1 F t (2) 



N\- 



-W0-, S)-S(2, 3))G!(1, 2)1, 



where the last term proportional to the single-particle Green's function 
originates from the time derivative of the time-ordering operator. Requiring the 
continuity equation for the density and current deviation operators to hold leads 
eventually to the result (7.9). 

The generalized Ward identity (7.9) can also be rewritten in terms of the 
scalar and vector vertex functions through eqs. (C.8) and (C.9). Projection from 
the left and the right onto a pair of Gf 1 with suitable arguments yields in fact 

(7.11) A r(1; 2; 3) + v 3 • T(l, 2; 3) = Ui(2, 3) - 3)) Gt\1, 2) = 
au n 



= j(S(2, 3)-£(l, 3)) 



V 3 J. h * T72 

3<! 2m 



2)-M(l, 2) 



(cf. eqs. (B.4) with U=0, (2.2), (3.11), and (3.12)). Notice that the Hartree 
contribution to the self-energy operator (as well as the potential term) drops out 
from the right-hand side of eq. (7.11). This remark suggests that eq. (7.11) 
actually holds for the irreducible vertex functions introduced in appendix C. We 
can, in fact, express 

(7.12a) 7X1, 2; 3) = f(l, 2; 3) + |d45f(l, 2; 4)w(4, 5) z (5, 3) , 

(7.126) T(l, 2; 3) = f(l, 2; 3) + Jd45f(l, 2; 4)v(4, 5) *(5, 3) , 
that give 

(7.13) J-r(l, 2; 3) + V 3 • 2; 3) = -^-f\l, 2; 3) + V 3 • f(l, 2; 3) + 



+ 



| d45r(l, 2; 4)v(4, 5) 



dt. 



■ z (5, 3) + V 3 -*(5, 3) 
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The continuity equation (6.9) for the induced density and current within linear 
response together with the symmetry properties of the correlation functions 
under the interchange of their arguments, guarantee the expression within 
brackets to vanish and allow us to replace the vertex functions in eq. (7.11) by 
their irreducible counterparts. 

The ordinary Ward identities, which express a constraint between the 
derivatives of the self-energy operator and the vertex functions in the limit of 
long wave lenghts and small frequencies [9], can now be recovered from the 
generalized Ward identity (7.11) as follows. Let us consider, in particular, a 
translationally invariant system that admits the Fourier representations (D.13). 
After elimination of the spin variables eqs. (7.11) and (7.13) reduce thus to 

(7. 14) q l\p; q)-q-f(p;q) = q -^q-p + ±- [M(p - q/2) - M(p + q/2)] , 

m n 

where p = (p, p ) and q = {q, q ). Taking the limit as q-+ (for given p) of eq. 
(7.14) yields 

(7.15a) r(„;0) = l 3M(P) 



3hpo 



when q/q ^>Q, and 



(7-156) f(p; 0) = ^ + V* p M(p), 

m 

when g/g ->oo. Notice that, on the Fermi surface, the Ward identities (7.15) 
relate the renormalized vertices to the wave function renormalization 



(7.16) 



1- 



BM(p) 
dhpo 



and to the effective mass of a quasi-particle 



(7.17) 



m 



1 + 



m 



h\p T \ dh\p\ 



M(\p\, sf) 



Equation (7.15a), in particular, is needed, e.g., to derive Landau's transport 
equation for a Fermi liquid from the Bethe-Salpeter equation [9]. 

Similar identities have been utilized in the theory of the inhomogeneous 
interacting electron gas to reveal the short-range density dependence of the 
mass operator [17], as well as in the theory of shallow impurities states in 



APPLICATION OF THE GREEN'S FUNCTIONS METHOD ETC. 



25 



semiconductors [18]. For these inhomogeneous systems the identity 
corresponding to eq. (7.15a) can be obtained alternatively in the following 
way [19]. We may regard the scalar potential U in eq. (3.14) as a slowly varying 
perturbation in time which is uniform throughout the system and set accordingly 
£7(1) = -aq with a infinitesimal. This perturbation induces a change in the 
inverse single-particle Green's function which can be written in two ways. On the 
one hand, we can make use of the scalar (irreducible) vertex function (3.17) to 
write 

(7. 18) SGT Hi, 2) = - \ d3TU, 2; 3) 8V(3) , 
with (cf. eq. (3.15)) 

(7. 19) 57(3) = -aq-ihj d2v(3, 2) 2 + ) . 

On the other hand, we can consider the shift (7.19) as generated by the gauge 
transformation (7.1) with A(l) = act u c being the velocity of light, and write 
according to eq. (7.6) 

(7.20) 8Gi(l, 2) = iqa(t 1 - t 2 )G 1 (l, 2) . 

n 

Equation (7.20) yields (cf. eqs. (A. 7), (B.2), and (B.4)) 



(7.21) SGT\1, 2) = p*(t 1 -WGi 1 (.l, 2) = q* 



2)-l(t 1 -t 2 )M(l, 2) 
n 



as well as 87(3) = - «q, since SGi(2, 2 + ) = 0. Comparison of eq. (7.21) with eq. 
(7.18) leads to the Ward identity 

(7.22) Ut, - t 2 )M(l, 2) = 2) - / d3f(l, 2; 3) . 



Upon taking the time Fourier transform, eq. (7.22) reads 

(7.23) gfo r 2! Q ) = ^ ~ r 2) - / dr 3 An, r 2 ; r,\Q, *>-> 0) , 

which reduces to eq. (7.15a) for a homogeneous system. 

From eq. (7.9) it can also be shown that any «$-derivable» approximation 
satisfies an important sum rule known as the Thomas-Reiche-Kuhn sum rule (cf. 
appendix E). 
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8. - Optical properties of semiconductors. 

The purpose of this section is to apply the formal tools developed so far to the 
discussion of the optical properties of semiconductors. In particular, we shall 
combine the treatment of the dynamics of the many-body effects beyond the 
one-electron model with the geometry characterizing a crystalline 
semiconductor. After some general considerations which apply to 
semiconductors (or insulators) with arbitrary point-group symmetry and which 
can be formally carried out to all orders of (many-body) perturbation theory, we 
shall present numerical calculations for (cubic) covalent materials whereby the 
relevant many-body effect is the excitonic effect. 

8'1. Maxwell's equations and local-field effects. - Charge densities and 
currents (over and above the ground-state values) can be partitioned into 
externally specified and induced contributions. The corresponding Maxwell's 
equations then read 

V-E(r, t) = 4-( ?en (r, t) + PM (r, t)) , 

V B(r, = 0, 

(8.1) \ 3 

VxE(r, t)=-fB(r, t) , 

at 

V x B(r, t) =^(J e Jr, t) + J ind (r, t)) + ± J-£(r, t) . 

c c at 

Here the induced charge density and current can be obtained within linear 
response by multiplying the induced part of the number density (5.4) and current 
(5.5) by the electronic charge (q = -e). A comment about the nature of the scalar 
and vector potentials that enter these expressions is in order at this point. 
Although the choice of the gauge has been left unspecified in sect. 5 as well as in 
appendix C, the forms (2.1) and (C.l) for the unperturbed and interaction parts 
of the Hamiltonian, respectively, are actually consistent with the Coulomb (or 
transverse) gauge, in which 

(8.2) V-A(r, <) = 0. 

The scalar and vector potentials entering eqs. (5.4) and (5.5) are thus to be 
understood as the external scalar potential and the total (transverse) vector 
potential which includes the induced contribution [20]. 

The Coulomb gauge yields a natural separation between the transverse (T) 
and longitudinal (L) parts of the electric field. We find it also convenient to 
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introduce the so-called perturbing electric field [21] 



(8.3) E p (r, = - - ^-A(r, t) - V P ext (r, t) = E ext (r, t) + Ef%r, t) = 



= E(r, t)-Er(r, t) , 



in terms of which the induced charge density and current can be expressed. To 
this end, we need to exploit the generalized Ward identity (7.9) and the 
symmetry condition (6.3), to obtain 



(8.4) 



and 



(8.5) 



■ i^X^Tu r 2 ; w) + Vj • z E (r b r 2 ; w) = , 

i*XB.(ru r z> w) + Vi • zbGi, r 2 ; co) + — (N\p(r 2 )\N) Vi £(rj - r 2 ) = , 



m 



toxR(r!, r 2 ; w) + V 2 • ^(^i, r 2 ; <u) = , 



^Ze(»-i, r 2 ; o>) + V 2 • ^(n, r 2 ; w) + — (N\p(rJ\N) V 2< ?(r 2 - = . 



In eqs. (8.4) and (8.5) advantage has been taken of the time translational 
invariance of the unperturbed system by introducing the time Fourier 
transform, and the analytic continuation from time-ordered to retarded 
correlation functions has been performed according to the prescriptions of sect. 
5. Recall also that eqs. (8.4) and (8.5) hold for any conserving (^-derivable) 
approximation to the correlation functions. 

Equations (8.4) ensure that the induced charge density and current obtained 
from eqs. (5.4) and (5.5) satisfy the continuity equation. Equations (8.5), on the 
other hand, enable us to rewrite 



(8.6) 



PincA 



(r; w ) = - %- f dr' * R (r, r'; co) • E F (r'; co) , 

Ico J 

JmdT', a) = - iw j dr' a(r, r'; co) ■ E F (r'; <o) , 
where we have introduced the so-called quasi-polarizability tensor 



(8.7) «(r, r'; *>) = - ^ (xr(t, r'; co) + ± (N\p(r)\N) S(r - r') 1 j , 

1 being the unit dyadic. In deriving eqs. (8.6) the surface contributions that 
originate from an integration by parts have been assumed to vanish. 
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It is common practice when dealing with dielectric materials to implement 
Maxwell's equations by introducing the polarization and displacement vectors 



(8.8) 



P(r; to)= - —Jind(r; co) . 

ICO 



D(r ; co) = E(r; co) + AnP(r; co) . 



(We assume throughout the absence of magnetic phenomena.) Notice that in eq. 
(8.8) E(r; co) is the total electric field, while in eqs. (8.6) only the perturbing 
electric field (8.3) enters, thereby justifying calling the tensor (8.7) the quasi- 
polarizability. 

The Coulomb gauge allows us also to express the scalar potential in terms of 
the instantaneous charge density and the vector potential in terms of the sole 
transverse currents, as they satisfy the equations [22] 



(8.9) 



V 2 ?(r; co) = - 4^( Pext (r; co) + PM (r; co)) , 



V 2 + ^jU(r; co) = -^(Jf xt (r; co) + JfV; co)) . 



In what follows we shall assume that Jf xt = inside the material. 

All equations considered so far hold for the microscopic fields which exhibit 
large and irregular variations on the atomic scale. This phenomenon occurs in a 
crystal irrespective of the fact that the external fields may vary only on a 
macroscopic scale corresponding to several crystal cells. These features (which 
are called local-field effects) thus render the response to an electromagnetic field 
of a crystal different from that of a homogeneous material. The difference can be 
expressed by stating that any response function of the crystal satisfies the 
(geometrical) property 



(8.10) 



f(r + n, r' + n; co) =/(r, r'; to) 



for any vector n of the Bravais lattice; for a homogeneous system, on the other 
hand, there is no restriction on the translation that appears at the left-hand side 
of eq. (8.10). 

To fully exploit the symmetry property (8.10), it is convenient to introduce 
the space Fourier transform, which reads 



(8.11) f{r,r';to) = 



1 



£ 2 2 exp [i(k + G) ■ r]f(k + G,k + G'\ co) exp [- i(k + G') • r'] . 

" * GC' 
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Here Q is the volume occupied by the crystal, A: is a Bloch wave vector confined 
to the first Brillouin zone (BZ), and G and G' are reciprocal lattice vectors. The 
appearence of the same k in the two exponential factors at the right-hand side of 
eq. (8.11) enables us to consider k itself, besides w, as a parameter in the Fourier 
transform of eqs. (8.6)-(8.9). 

Microscopic fields are not, however, the quantities which are dealt with in 
ordinary electrodynamics where one considers field quantities that vary on a 
macroscopic scale and are thus experimentally accessible. To obtain the 
macroscopic from the microscopic fields an averaging procedure is necessary 
which has the result of smoothing out the irregular fluctuations of the 
microscopic quantities [23]. A suitable averaging procedure for functions which, 
once represented by the Fourier series 

1 BZ 

(8. 12) <7(r; w ) = -L_ ^ S exp [i(k + G) • r] g(k + G» , 

contain only small Bloch wave vectors (i.e., \k\ « \G\), is to take the average over 
a unit cell of the cell-periodic part only [24] 

1 f 1 BZ 

(8. 13) g(n; = — drgr(r; «) = — 2 exp [ik ■ n] g(k + 0; to) . 

Uc w y/o k 

Here Q c [n] is the unit cell centred at n, Q c is its volume, and the approximation in 
the last line neglects the variation of exp [ik • r] over the unit cell. Consistently, 
one may replace n in eq. (8.13) by the continuous variable r. 

The averaging procedure (8.13) suffices to eliminate the rapidly varying 
fields from eqs. (8.8) and (8.9), provided the external fields are slowly varying 
over the unit cell, but it fails, in general, for eqs. (8.6) where the average of a 
product needs to be taken. In practice, however, we can neglect the microscopic 
components {G 0) of the perturbing electric field for values of the frequency u> 
restricted to the optical range [21]. In fact, combining the Fourier transforms of 
eqs. (8.8) and (8.9) yields 

(8. 14) E T (k + G; oj) = -I—LDtOc + 

(k + GY 

which tipically gives (Cr#0) 

/oic\ \ E T(k + G; <o)\ 

\E T (k + 0; to) I 

for h(o~5 eV. Similarly, the microscopic components of the longitudinal part of 
the displacement vector can be neglected, since one may show that 

(8. 16) D L (k + G; w ) = Ei*(k + G; w ) . 
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8'2. Macroscopic dielectric tensor. - The optical properties of semiconduc- 
tors can be conveniently described in terms of the macroscopic dielectric tensor 
e M (k; &>) which relates the macroscopic components of the displacement vector 
and of the total electric field: 

(8. 17) D(k + 0; to) = e M (k; <u) • E(k + 0; w) . 

e M (k; oj) can, in turn, be expressed in terms of many-body response functions by 
relating it to the macroscopic component a(& + 0, A; + 0; &>) of the quasi-polariz- 
ability tensor (8.7). To this end, we need first to relate the macroscopic 
components of the total and of the perturbing electric fields via 

(8. 18) E(k + 0; w ) = [T - A-(kk) ■ a(k + 0, k + 0; w )] • E*(k + 0; w ) , 

where (kk) is the longitudinal projection operator. Inverting this equation in 
favour of E F (k + 0; w) and using eq. (8.8) we find eventually 

(8.19) e M (A:; w) = "T + 4?za(k + 0, k + 0; <*>) -ft - A7z{kk) ■ *(k + 0, A; + 0; w)]" 1 . 

The matrix within brackets at the right-hand side of eq. (8.19) may then be 
explicitly inverted, to give [21] 

(8.20) e M (k; w) = T + 4-a(Ar + 0, k + 0; co) ■ 

| Mkk) • a(fc + 0, A: + 0; co) 
1-iTzk- a(k + 0, k + 0; a>) • k_ ' 

There still remains to evaluate the left and right longitudinal contractions of 
the tensor a(Ar + 0, k + 0; w). This can be readily done by taking the space 
Fourier transform of eqs. (8.5b), (8Ab), and (8.4a). We obtain, in the order, 

Q 2 

(8.21a) a(A? + 0, A; + 0; w) ■ k = - -V zr(* + °, k + 0;co), 

oik 

* „ a 2 
(8.216) k ■ *(k + 0, k + 0; co) = - \ x^k + 0, k + 0; w ) , 

a 2 

(8.21c) • a(Ar + 0, A; + 0; &>) • = - — ynik + 0, A: + 0; a>) . 

Insertion of eqs. (8.21) into eq. (8.20) allows us to calculate the macroscopic 
dielectric tensor by many-body techniques. 

In eq. (8.20) the parameter A; is assumed to be restricted in such a way that 
|A;| _1 is much larger than the characteristic microscopic lenght scale, i.e. the 
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dimension of the crystal cell. The frequency co, on the other hand, ranges over 
the optical region and thus it is quite comparable with the semiconductor band 
gap. One may then envisage eliminating k from eq. (8.20) by formally taking the 
limit as k^>0, while keeping « finite. In performing the limit, however, some 
care should be exerted owing to the divergence of the (Fourier transform of the) 
bare Coulomb potential in this limit. We shall show that the macroscopic 
dielectric tensor is well-behaved in this limit, irrespective of the divergence. 

To single out from the correlation functions the part containing the 
divergence of the bare Coulomb potential, it is convenient to rearrange the 
integral equations they satisfy (cf. eqs. (3.24) and (C.31) for the time-ordered 
counterparts) in the following way. We first introduce the proper correlation 
functions (denoted by a circumflex) as being the solutions to the equations 

(8.22a) xn(k + G,k + G'; «) = ■ 
= X E (& + G, k + G'; co) + 2 X u(k + G,k + G"; co)v(k + G") X v.(k + G", k + G'; co) , 

G"#0 

(8.226) Mk + G,k + G'; co) = 

= %(k + G, k + G'; co) + 2 *r(* + G,k + G"; co) v{k + G") j R (k + G", k + G';co), 

(8.22c) £t(* + G,k + G'\ co) = 
= %(k + G, k + G'; w)+2%(k + G,k + G"; co)v(k + G")xn(k + G", k + G'; co) , 

(8.22d) %{k + G,k + G'; co) = %(k + G,k + G'; co) + 

+ E xn(k + G,k + G"; co) v(k + G") j R (k + G", k + G'; co) + 

G"#0 

+ 2 2%(k + G,k + G"; co)v(k + G") ■ 

• xdk + G", k + G'"; co) v(k + G'") %{k + G'", k + G'; co). 

Here v(k + G) = An#l\k + Gf and the retarded irreducible correlation functions 
(denoted by a tilde) are obtained by analytic continuation from the time-ordered 
counterparts introduced in appendix C. Owing to the absence in eqs. (8.22) of the 
factor v(k + 0) which would be the source of the divergence in the limit 0, the 
proper correlation functions can be interpreted as the correlation functions of a 
hypothetical medium whereby the interparticle potential is replaced by a 
Coulomb potential with a small-momentum cutoff [25]. These functions are thus 
well-behaved in the limit »0. 

The G = G' = elements of the correlation functions that enter the 
expression (8.20) of the macroscopic dielectric tensor can be cast in a form which 
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isolates the divergent factor v(k + 0) explicitly [26]: 

Xr(A + 0, k + 0; w ) 



(8.23a) XR (k + 0,k + 0;co) = ~ „ 

1 - v(Ar + 0)xr(At + 0, A + 0; &j) 

(8.236) Jr(* + 0, h + 0; a.) = 

1 - v(k + 0)xr(k + 0, k + 0; &>) 
(8.23c) x*(k + 0, /c + 0; &j) = 

v(ft + 0)-f E (fc + 0, fc + 0; w) 



1 — v(k + 0) xe(k + U, ft + 0; to) 

(8.23d) x R (ft + 0, k + 0; w) = 

C c v(fc + 0)^ + 0, k + 0; w) 

1 - v(ft + 0)xR\k + 0, ft + 0; to) 

We are now in a position to rewrite the macroscopic dielectric tensor (8.20) in 
a form suitable for taking the limit A; ->0. Inserting eqs. (8.23) into eqs. (8.21) 
and into the Fourier transform of eq. (8.7) and entering the results into eq. 
(8.20), we obtain after straightforward algebraic manipulation: 

(8.24) e M (fc; w ) = ( 1 - ^) 1 -^f^(ft + 0, k + 0; w ) , 

which is indeed well-behaved in the limit ft— > 0. Equation (8.24) is formally valid 
to all orders of (many-body) perturbation theory and for crystals of arbitrary 
point-group symmetry [27]. For cubic materials eq. (8.24) takes a particularly 
simple form, to be discussed next. 

8'3. Cubic materials. - For materials with cubic symmetry the macroscopic 
dielectric tensor (8.24) is equivalent to a scalar in the limit k^>0, i.e. it is 
diagonal with equal elements: 

(8.25) lime M (ft; eo) = e M (w) 1 . 

A:—. 

Here e M (w) can be obtained, in particular, from the longitudinal-longitudinal 
element of the matrix (8.20) by taking into account eqs. (8.21c) and (8.23a): 

(8.26) e M (w) = 1 - Hm v(k + 0) X n(k + 0, k + 0; w) , 
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where 



(8.27) 



Xz(k + 0, k + 0; w) fc^b k 2 XR(co) , 



for e M (A;; to) to be well-behaved in this limit. 

Equations (8.25) and (8.27) can be proved by exploiting the Lehmann 
representation (5.10) of the proper correlation functions (8.22d) and (8.22a), in 
the order. For completeness, we sketch the proof below. 

i) Proof ofeq. (8.25). Equation (8.25) results from the following symmetry 
property of the proper current-current correlation function: 

(8.28) m,' 1 ■ %{&k, fflk; *>)•&,= %{k, k; to) , 

where &l is the rotation part of a symmetry transformation {0l\w} of the 
crystal. Equation (8.28), in turn, follows from the invariance under this 
transformation of the set of projection operators (built from many-body 
eigenstates of given energy) that enter the Lehmann representation. Taking the 
limit as of both sides of eq. (8.28) and choosing successively ffi to be, e.g., a 
rotation by tz about the a?-axis, a rotation by n about the y-axis, and a rotation by 

2;r/3 about the Ill-axis [28], one may readily prove that xr(0, 0; <u) is equivalent 
to a scalar. 

ii) Proof of eq. (8.27). Equation (8.27) follows from the Lehmann 
representation of the proper density-density correlation function together with 
the property 



for s¥=0 (cf. eqs. (5.10) and (5.11)) and the fact that the system possesses an 
energy gap. (QED) 

The result (8.25) simplifies considerably the discussion of the optical 
properties of cubic semiconductors. In particular, one may write for the induced 
current (cf. eqs. (8.66), (8.16), (8.18), and (8.20)) 



where the transverse and longitudinal components have been explicitly indi- 
cated. Equation (8.30) implies, in particular, the well-known results that there is 
no transverse screening effect in the long wavelength limit and that a transverse 
perturbing electric field excites only transverse induced currents [25]. It then 



(8.29) 




(8.30) 
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follows that we can find purely transverse solutions to Maxwell's equations (8.9), 
whereby ?(r; w) = 0. Setting 

(8.31) A(r; u) = A exp[i/-r] 
and taking 

(8.32) J?*(r; w ) = - J ^L A{r . w) 



4;: / c — 

we find that the propagation vector A satisfies the dispersion relation 

(8.33) ^ = 4e M M. 

cr 

In writing eq. (8.32), we have assumed that the real and imaginary parts of A are 
small enough for the space Fourier transform of eq. (8.31) to contain only small 
wave vectors A;. In this way the same constant e M (w) is pertinent to all Fourier 
components of the wavepacket (8.31). 

Conventionally, one writes the propagation vector A in terms of the index of 
refraction (n) and of the extinction coefficients (x) 

(8.34) A = s-(n + ix), 

c 

where s is the unit vector in the direction of propagation. One can then express 
the average energy flux in the medium corresponding to the (real part of the) 
solution (8.31) through a Poynting's vector of the form 



(8.35) (S) = §-fn^A 2 exp 

07Z C 



o CO - 

2 — xs • r 



The coefficient of the exponentially decreasing factor in eq. (8.35) is called the 
absorption coefficient r, [39], and may be also expressed in terms of the 
imaginary part e 2 M of the macroscopic dielectric function e M (&j) by using the 
dispersion relation (8.33): 



(8.36) >;(«) = - 



oi e 2 (co) 



c n(e») 



This expression for -q can alternatively be obtained by calculating the power 
absorbed by the medium. To this end, we identify the rate at which the 
electromagnetic field does work on the system as the time rate of change of the 
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mean value of the system Hamiltonian (2.1) [24] 

(8.37) dW=A.( Y (t)\H\nt)) =i(nt)\[A'(t), (H + H'{t))]\Y(t)) , 

where \T(t)) is the wave function of the system driven by the electromagnetic 
field and H'(t) is the interaction Hamiltonian (C.l) (without the diamagnetic 
term). Standard manipulations lead to 



(8.38) ^= - { dr V xt (r, t) ■ J ind (r, t) + 

+ I/drA( M ).| 



JUr,t) + ^(N\p(r)\N)A(r, t) 



where use has been made of the continuity equation and a surface term 
originating from integration by parts has been dropped. To recover from eq. 
(8.38) the classical expression involving only the induced current and the 
(perturbing) electric field, we should perform the time average of eq. (8.38) over 
the time interval (- T, +T) during which the perturbation acts on the system. 
(For periodic disturbances, it is actually enough to average over one cycle.) We 
obtain 



■T 



where we have performed an integration by parts over t' and recalled eq. (8.3). 

Equation (8.39) holds quite generally for any type of material. For a cubic 
semiconductor we may consider a solution of the form 



(8.40) 



E p T (r; co) = rc[E p T (r) S(co - w ) + El(r)* $(«> + o> )] , 



EtW = i~^o exp \ik ■ r] , 



(ftj >0), for which 

(8.41) l^y^^\dr\E\{r)\\ 

Defining the absorption coefficient as the ratio between the (time-averaged) 
energy absorbed by the system (per unit time and per unit volume) and the 
energy flux in the medium (for the same unit volume), we are able to recover eq. 
(8.36) from eq. (8.41). 
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One may also verify that the (time-averaged) energy absorbed by the system 
per unit time can be obtained directly for a transverse vector potential 
corresponding to eq. (8.40) by calculating the transition probability per unit time 
from the Golden rule. This calculation shows, in addition, that 



for the transverse submatrices in the limit A;— » 0. This result can be interpreted 
by saying that the proper density-density correlation function describes the 
response of a (cubic) semiconductor (or insulator) to a transverse electric field of 
long wavelength [25]. Notice that eq. (8.42) could have alternatively been 
obtained from eq. (8.23d) by showing that the second term at the right-hand side 
of that equation vanishes by symmetry in the long-wavelength limit [30]. 

8*4. Time-dependent-screened-Hartree-Fock approximation in a local- 
orbital basis. - The result (8.26) for the macroscopic dielectric function holds 
quite generally whatever form of the proper density-density correlation function 
is inserted at its right-hand side. We now discuss some numerical results which 
have been obtained for (cubic) covalent crystals such as diamond [31] and 
silicon [32, 33], and more recently for GaP [34], whereby the important many- 
body mechanism is the screened interaction within an electron-hole pair created 
in the polarization process. Parallel calculations of the photoabsorption cross- 
section in atomic systems within the RPAE approximation have also emphasized 
the importance of the above mechanism [35], although methods for different 
systems rely on different computational schemes. In particular, for 
semiconductors the use of a local-orbital basis to represent single-particle wave 
functions has made it possible the inclusion of the excitonic effect on top of local- 
field effects. 

As discussed in sects. 2-4, the dynamics of the many-body processes 
occurring in the system response is embodied in the choice of the functional form 
of the self-energy operator, since the effective two-particle interaction can be 
obtained from it by functional differentiation. Consistency with conservation 
laws further restricts the choice of the self-energy, as discussed in sect. 6. For 
the purpose of describing the optical properties of semiconductors it is 
apparently sufficient to take for the non-Hartree part of the self-energy operator 
(cf. eqs. (3.11) and (3.12)) 



as well as to neglect the variation of the dynamically screened interaction while 
calculating the functional derivative: 



(8.42) 



teik + 0, A: + 0; W ) T = + 0, k + 0; W ) T , 



(8.43) 



Mi\, 2) = ihW(l + , 2)G 1 (1, 2), 



(8.44) 



SMiX, 2) 
8GA 4) 



= 3) £(2, 4)W(1 + , 2). 
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Equation (8.44) would become exact if W in eq. (8.43) would be replaced by the 
bare Coulomb potential. In this case the system response is said to be 
approximated within the time-dependent-Hartree-Fock approximation, which 
has been used with success for several atomic systems. For semiconductors, 
however, the screening of the electron-hole attraction is an essential feature 
which cannot be neglected in calculating the optical spectra. The same conclusion 
will consistently be reached when discussing the single-particle energies (sect. 9) 
and the (bound) excitonic levels (sect. 11). 

The form (8.43) for the mass operator (the so-called GW approximation) has 
been first discussed by Hedin in conjunction with properties of the electron 
gas [36]. This form is consistent with conservation criteria provided that the 
dynamically screened interaction is evaluated within the random-phase- 
approximation (RPA), i.e., by approximating the (irreducibile) scalar vertex 
function that enters the expression (3.26) of the irreducibile polarizability by its 
lowest-order expression. An example is given in fig. 3e). Going beyond the RPA 
approximation (by introducing, e.g., a screened electron-hole interaction in the 
bubbles) requires, in principle, the addition of other terms to the mass operator 
in order to be consistent with conservation criteria (cf. fig. 3/)). By the same 
token, the approximation (8.44) where some terms at the right-hand side have 
been dropped, violates in principle the conservation criteria. In practice, 
however, in calculating the optical response one usually approximates the 
dynamically screened interaction by a static screened potential obtained from a 
phenomenological dielectric function 

(8-45) W(l + , 2) - Wfru r 2 ) S(tt - W , 

and keeps the form (8.44) for the kernel of the Bethe-Salpeter equation. 
Consistency with conservation criteria is then limited to a numerical check on the 
accuracy to which the longitudinal sum rule is satisfied. Further comments on 
the practical importance of a formal violation of conservation criteria in the 
calculation of the band structure of covalent semiconductors will be given in 
sect. 9. 

We pass now to solve the integral equation (C.16) whose kernel (C.15) is 
approximated by eqs. (8.44) and (8.45). It is convenient to eliminate at the outset 
both the spin variables according to the prescriptions of appendix D and the time 
variables by first pairing them as in eq. (C.18) and then taking the Fourier 
transform with respect to the resulting relative time. We obtain 

(8.46) L(rj, r 2 ; r[, r' 2 \o>) = L (r u r 2 ; r[, r' 2 \o>) + 

+ ~j dr s dr 4 L (r u r 4 ; r[, r z \u>)W s {r z , r 4 )L(r 3 , r 2 ; r 4 , r' 2 \<o), 
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where we have introduced the notation 

(8.47) L (r u r 2 ; r[, r^) = 2 j —G^r,, r£ a> + a/) G,(r 2 , rj; a/) 

for the noninteracting part. Notice that the factor of two at the right-hand side of 
eq. (8.47) originates from the spin degeneracy. Notice also that in eq. (8.46) the 
frequency to is a parameter. 

To proceed further, we represent all single-particle Green's functions 
entering eq. (8.46) in the quasi-particle approximation of the form 

(8.48) Gjfo, r 2 ; w)~2 



k ha,- s k + iosgn (sfc - e F ) ' 



where <?— >■ + and the Fermi level s F is placed within the gap. In eq. (8.48) the 
u k (r) are (orthonormalized) single-particle wave functions which are assumed to 
be known by a previous solution of a band eigenvalue problem. In practice, they 
are approximated by a or a local-density calculation. Although eq. (8.48) is not 
the most general form for a single-particle Green's function, it turns out to be a 
suitable approximation to describe the optical properties of semiconductors. In 
sect. 9 we shall justify the choice (8.48) as well as its internal consistency. 
Entering eq. (8.48) into eq. (8.47) yields 



(8.49a) L (r u r 2 ; r[, ro\u) ■ 



2 u kl (r{) u^Loikukz; k[, k' 2 \u>) u ki (r[)* u ki (rQ* 



where 

(8.496) L Q {k u k 2 ; k[, k 2 \oj) = 



d(e k „ - e F ) d(e F - e k 6(e. k > 2 - s F ) 6(e F - SjJ 



s*i - £*„ + + htn £*- 2 - e k2 - iS + ha 



Here the step functions discriminate conduction bands from valence bands 
states. (Notice that the infinitesimal imaginary parts in the energy denominators 
in eq. (8.496) correspond still to time-ordered quantities. Analytic continuation 
to retarded quantities will eventually be performed by replacing i<?— > — iS in the 
second energy denominator within brackets). We may similarly expand 

(8.50) L(r b r 2 ; r[, r 2 \co) = £ i6 tl (ri)^ 2 (r 2 )L(fc 1 ,/c 2 ; k[, k^) u^ird* u^)* , 

il&2 k'i k% 



and transform eq. (8.46) into an integral equation for the coefficients at the right- 
hand side of eq. (8.50) by invoking the assumed orthonormality of the single- 
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particle wave functions. A more handy matrix equation can instead be obtained 
by representing the single-particle wave functions in eqs. (8.49a) and (8.50) in a 
local-orbital basis [31]: 

(8.51) u nk (r) = " 1/2 2 exp [ik ■ I] 0„(r - 1 ) c„ K (&) . 

Here we have specified the label k to signify the pair (n, k) where n is a band 
index and A; is a Bloch wave vector, Jf is the number of lattice sites, and the 
orbital £„(r- i) is localized about the lattice site I. In eq. (8.51) it is understood 
that the set is sufficient to describe the valence and conduction bands which 
are coupled by optical excitations. 

The following steps allow us then to transform the Bethe-Salpeter equation 
(8.46) into a matrix equation: 

i) set r x = r[ and r 2 = f-2 in the light of eq. (C.18); 

ii) multiply both sides of eq. (8.46) by Q' 1 exp [- i(g + G) • r 2 ] • 
• exp [i(q + G')-r 1 ], where q now stands for a Bloch wave vector, and integrate 
over r x and r 2 ; 

iii) use the symmetry properties of the screened interaction W s (r 3 , r 4 ) 
when shifting both coordinates by the same lattice vector as well as when 
interchanging r 3 «±r 4 . 

After suitable shift of integration and of (discrete) summation variables, we 
obtain eventually 

(8.52) -f 2 A c (q + G) \n cc .(q; *>) - n°Aq; *») - 

"0 c,c' [ 

- f 2 ^cAq; *») V?Aq)N°rAq; ")U C % + G') = 0. 

In eq. (8.52) Q is the volume of the Wigner-Seitz cell, the collective index «c» 
stands for the triad (I, v, p) which appears in the form factors 

(8.53) A c (q + G) = \ dr #*r) exp [- i(q + G)-r] &(r - /) , 
and we have introduced the matrices 

1 BZ 

(8.54) NAq; = 2 2 c* 2 (k)c^(k + q)exp[i(k + q)-l]- 

■L(nsk', n^ + q; nje' + q, nzk\u) exp[- i(k' + q) • l']cp ni (k' + q)c^(k') , 
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I BZ 

(8.55) NUq; «) = 2 2 + g) exp + «)•(/- /')] • 

^ ' k 



2z 
ft 



e( £?t2 (fc) - £ F ) 8( £F - + g)) 
s n Xk + q) - tjk) + iS - hoj 



d( SF - Sn2 (k))e( Sni (k + q)-e F y 



Cs4k)cf ni {k + q), 



e Ml (Ar + q) - e^k) - iS - hco 
(8.56) Vg.(g) = 2)exp[-tg-m]- 

m 

• <tf(r' - m) &(r - 1 - m)\W s (r, r')| (r') &.(r -!')>, 

which represent the irreducible polarizability, its RPA approximation, and the 
screened interaction in the local-orbital basis, respectively. For reasons that will 
soon become clear, eq. (8.56) is called the (screened) exchange Coulomb matrix. 

Equation (8.52) can, in general, be satisfied if the matrix within braces 
vanishes identically. This is the desired matrix equation, which can be formally 
solved to give 



(8.57) 



N = N° 



2 



in matrix notation. The practical convenience of this computational scheme is 
apparent at this point whenever the spatial localization of the orbitals <f> 
drastically limits the range of the index c in eq. (8.52) (cf. eq. (8.53)). In the 
limiting case of extreme localization one may, in particular, recover the classical 
Clausius-Mossotti relation from eq. (8.57) and below (cf. appendix F). 

Within the above approximations, the irreducible polarizability takes on the 
form 

(8.58) x {q + G, q + G', w ) = -±- £ A c (q + G)(- ih)N cc (q; u)A*.(q + G') . 

"0 cc' 

The proper polarizability, which is needed to evaluate the macroscopic dielectric 
function (8.26), may be similarly expressed as 

(8.59) x(q + G,q + G';<o) = ±- £ A c (q + G)S c Aq; «>)A*(q + G') , 

"0 cc' 

where S is called the (proper) screening matrix. S can be determined by entering 
eqs. (8.58) and (8.59) into the integral equation (8.22a). One finds 



(8.60) S = (- ih)N[l + ihV c NY l = (- ih)N° 



! + (f V s *- V c )(-ih)N° 
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in matrix notation, where in the last line use has been made of eq. (8.57), and 
(8.61) V c Aq) = 2 A*(q + G) v(q + G)A,(q + G) 

"0 G*0 

is the (proper) Coulomb matrix. By adding to the right-hand side of eq. (8.61) the 
missing term with G = 0, one obtains the full Coulomb matrix which will be 
needed in sect. 9. The latter can also be cast in the form (8.56) with the 
interchange r+±r' between the first and the last orbital. This justifies calling eq. 
(8.56) the (screened) exchange counterpart of the Coulomb matrix. 

The limit as can be taken in eq. (8.59), whereby the property (8.27) 
results from eq. (8.53) and the orthogonality of orbitals referring to valence and 
conduction bands. 

Details of the calculations are given in refs. [31-34]. Here we only comment 
on the main features of the results. One is interested in comparing the 
experimental spectrum of e 2 (w) with the calculated spectrum obtained within 
alternative approximations based on the result (8.60), namely, i) the time- 
dependent-screened-Hartree-Fock approximation (TDSHF) that takes eq. (8.60) 
as it stands, ii) the random-phase-approximation (RPA) that neglects the 
excitonic effect and thus sets V sx = in eq. (8.60), and iii) the RPA with the 
further neglect of local-field effects (RPA) that sets also V c = in eq. (8.60) [37]. 
Figures 4 and 5 show the results for diamond and silicon, respectively. Common 
trends for both materials are that passing from the RPA to the RPA 
approximation furthers the discrepancy with experiment, while good agreement 
with experiment is re-established by the inclusion of the excitonic effect 
(TDSHF). Notice how these effects are more pronounced for diamond than for 
silicon, as expected from the different degree of localization of the orbitals. 
Similar results are found for GaP [34]. 




15 20 25 

energy (eV) 

Fig. 4. - Imaginary part of the dielectric function vs. energy for diamond (adapted from 
ref. [31]). The experimental spect rum is compared with three different calculations (see 
text). experiment, • • • RPA, <=oo RPA, x x x TDSHF. 
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5.0 3.0 
energy (eV) 

Fig. 5. - Imaginary part of the dielectric function versus energy for silicon (from ref. 
[35]). a) The experimental spectrum (■ • •) is compared with the TDSHF calculati on (qqq) ; 
b) comparison of three different calculations (see text). • • • TDSHF, oooRPA, aaa RPA . 



As a «rule of thumb», the excitonic effect tends to pull the strength of the 
absorption coefficient toward the low-energy side [38], a feature which can also 
be accounted for by the use of a simplified model [39]. The numerical results also 
show the appropriateness of describing the absorption above the fundamental 
gap in terms of localized orbitals, since large regions of A-space contribute to the 
main peaks. This has to be contrasted with the description of excitonic bound 
states below the fundamental gap whereby the use of the so-called effective- 
mass-approximation is found to be appropriate for semiconductors (cf. sect. 11). 



9. - Single-particle energy levels. 

In the previous section the underlying single-particle band structure needed 
for calculating the optical spectra of semiconductors (cf. eqs. (8.48) and (8.51)) 
has been assumed to be known and emphasis has been placed on two-particle 
effects over and above band-structure effects. Many-body techniques, however, 
allow us also to determine single-particle band structures of semiconductors. 
This study, on the one hand, is necessary for a consistency check on the optical 
spectra, and, on the other hand, has some interest of its own because single- 
particle energy levels can be directly probed whenever two-particle effects can 
be disregarded. In particular, assessing the role of many-body effects on the 
band structure of semiconductors bears on the so-called « band-gap problem » 
tiiat arises from the inability of conventional local-density calculations to 
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reproduce the experimental band gaps of semiconductors [40]. It turns out that 
the ansatz (8.43) for the mass operator is sufficient for an accurate description of 
the band structure of (covalent) semiconductors and wide-gap insulators, at least 
for the cases studied thus far [41-44], this finding being consistent with the 
relevance of the exeitonic effect on the optical spectra. 

The starting point for the connection between the single-particle Green's 
function and the band structure is the observation that the poles of the time 
Fourier transform of the single-particle Green's function (in the absence of an 
external time-dependent potential) correspond to the energies of excitation of 
quasi-electrons and quasi-holes [45]. As the poles might be complex, these are 
approximate excited states of the (N + l)-particle and (N - l)-particle system, 
respectively. A suitable algorithm to locate these poles starts by taking the time 
Fourier transform of Dyson's eq. (3.14) with [7 = 

(9.1) [hco - Kr)]Gi(r, r'; w ) - J dr\£(r, r"; co)G 1 (r J ', r'; w ) = S(r- r') , 

and then proceeds by considering the associated right and left eigenvalue 
problems [17] 

(9.2a) h{r) <p„{r; co) + J dr' 2 (r, r'; co) <pir'; co) = <£„(«) <p,{r; co) , 

(9.26) h{r) Ur; <o) + j dr' ft(r'; u)Z (r', r; co) = <?„(«) &(r; co) , 

whereby co is regarded as a parameter. In these equations, the self-energy acts 
as a nonlocal and energy-dependent potential [46]. The eigensolutions of eqs. 

(9.2) are assumed to form a complete set and to satisfy the orthonormality 
relation 

(9.3) J dr &(r ; «■>) Mr; «) = • 

The single-particle Green's function can thus be expressed in the bilinear form 

<p,fj; oj) <p„(r'; w) 



(9.4) G 1 (r,r';<o) = J J 



hco — <SSs>>) 



Equation (9.4) can be considerably simplified when the resonances occurring at 
hco ~ £„(<») are narrow. In this case, one usually adopts the so-called quasi- 
particle approximation [47] which assumes that the singularities originating 
from the energy denominator in eq. (9.4) form the dominant contribution to the 
single-particle Green's function: 

, Q „ r , . , v Ur; ReK})^(r'; Re{ M ,}) 
(9.5) dir, r';oj)=^X g, , 

no) — nco. 
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where &> v are determined by the condition 



(9.6) 



<?„( w v) = ha>,, 



and we have introduced the renormalization factor 



Since the skew-Hermitian part of the self-energy operator changes from being a 
positive operator to being a negative operator as hoi increases through the 
chemical potential [47] (which for a semiconductor lies within the energy gap s g ), 
the imaginary part of <?„ changes sign correspondingly. This remark implies, in 
particular, that if the self-energy did not depend on w, its skew-Hermitian part 
would vanish and consequently eq. (9.5) would become exact with g,, = l. This is 
indeed the case for the Hartree-Fock approximation 



where the factor of two in the Hartree term originates from the spin sum and the 
convergence factor in the frequency integrals selects the contribution of the 
states below the chemical potential in eq. (9.5). 

The Hartree-Fock approximation represents conceptually a good reference 
level to study many-body effects since it corresponds to the «best» single- 
particle picture whereby each particle is assumed to move under the influence of 
the average interaction with all other particles [48]. However, when applied to 
real materials it gives, e.g., energy gaps that are by far too large compared with 
the experimental values. For instance, in the case of diamond one obtains a value 
of 15.0 eV for the gap [49] which is about twice the experimental value (7.3 eV). 
The reason for this serious deficiency of the Hartree-Fock approximation can be 
traced in the lack of relaxation processes which should accompany the primary 
single-particle excitation, namely, the rearrangement of the passive electrons to 
the addition of an extra electron or to the creation of a hole. The success in 
describing the optical properties (sect. 8) suggests that the ansatz (8.43), which 
replaces the bare Coulomb interaction v of the Fock term with the dynamically 
screened interaction W, can actually account for the above polarization processes 
in semiconductors. Moreover, to address important questions such as the 
influence of the energy dependence of the dielectric matrix on the self-energy 
shifts and the goodness of the quasi-particle approximation in semiconductors, 



(9.8) Z HF (r, r') = - ih$(r-r')2 j dr"v(r-r") j 



^exp[te/<y]Gi(r w , r"; w ') + 



,71 
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we will keep the full time dependence of eq. (8.43) in the calculation and resort to 
a static approximation of the type (8.45) only for comparison. 

The adding of correlation effects within the GW approximation (8.43) on top 
of a Hartree-Fock band structure was the line followed in the original calculation 
for diamond [41] which was chosen as a prototype of covalent materials. The 
reasons for this choice were essentially the availability of accurate Hartree-Fock 
calculations for this light-core material as well as the appropriateness of a local- 
orbital representation for the valence and conduction bands which is needed for 
an accurate description of the screening properties (cf. subsect. 8'4). Later 
work [42-44] has used as starting point local-density band-structure calculations 
which have the advantage of being more readly available than Hartree-Fock 
calculations. Here we follow the approach of the prototype calculation for 
diamond [41] for pedagogical reasons. 

The local-orbital representation (8.51) for the single-particle wave functions 
allows us to reduce the quasi-particle equation (9.2) to an algebraic eigenvalue 
problem for the coefficients c^JJc) [50] 

(9.9) 2(<f*k*)k > + ilAMiJc; E n ik))W))c^k) = E n {k) c, n (k) . 

Here we have introduced the notation 

(9.10a) <jx|e(*)k> = 2exp[i*-l]<# ( ,(r)Mr)Mr-l)>, 

i 

(9.106) { l x\M(k;E)\ fJ i') = Zexv[ik-l](Mr)\M(r, r' ; E)\^(r' - I)) , 

i 

where h m (r) contains the Coulomb effects both of the nuclei and of the average 
electronic charge distribution in the ground state (Hartree term). Notice that 
the condition (9.6) has already been taken into account in eq. (9.9). 

For the expressions (9.10) to be of practical importance the summations over 
the lattice sites I should be rapidly converging. This is certainly the case for the 
matrix elements of the local Hamiltonian in eq. (9.10a), provided the local 
orbitals are sufficiently localized, while for the convergence of eq. (9.106) the 
short-range property in \r-r'\ of the mass operator in required in addition [17]). 
Crystal symmetry also considerably reduces the number of independent matrix 
elements in eqs. (9.10). For diamond a minimal set of ten independent matrix 
elements has proven sufficient [41]. 

The expression of the mass operator within the GW approximation to be 
inserted in eq. (9.106), namely 

+ CO 

(9.11) M(r, r'; E) = -±- \ dE' exvUE'flG^r, r'; E + E')W(r, r'; E') 

LTV J 
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can be split in terms of the Hartree-Fock M HF (the second term at the right-hand 
side of eq. (9.8)) and of a remainder M' , whereby the dynamically screened 
interaction in eq. (9.11) is replaced by 

(9.12) W'(r, r'; ^) = /dr 1 dr 2 v(r-r 1 ) ZT (ri, r 2 ; E)v(r 2 -r'), 

Xt being the (time-ordered) polarizability matrix. The matrix elements of M HF 
can thus be determined by fitting a Hartree-Fock band structure 
calculation [49]. The matrix elements of M', on the other hand, need to be 
explicitly evaluated. To this end, it is convenient to exploit the local-orbital 
representation for ^ T which is obtained from eq. (8.59) by replacing the proper 
screening matrix S with the full screening matrix S that contains also the term 
with G = 0. This is because, contrary to the description of the optical properties 
where the limit 9— »0 is in order, the short-range nature of the electronic 
correlations requires an accurate description of x T (<7 + G, q + G'; E) for all q. 
The result can then be expressed in terms of the matrix 

(9.13) W' m ^ m ;Aq; E) = 2 2 V m ,^(q)S^. v/ Aq; E)V v ., y , mW Aq) , 

where now the indices {a, z) refer either to a pair of bonding or antibonding 
orbitals, or to a bonding and an antibonding orbital. Details of the calculation are 

l r x 
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Fig. 6. - Quasi-particle band structure of diamond calculated within the Hartree-Fock 
approximation (dashed-dotted line) and within the GW approximation with the screening 
taken within the TDSHF approximation (full line) (from ref. [19]). 
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given in ref. [41] where the numerical methods adopted to make the scheme 
working in practice are also described. 

The quasi-particle valence and conduction band structure of diamond (that is, 
the real part of the eigenvalues E n (k) of eq. (9.9)) is shown in fig. 6 for two 
symmetry directions. The label TDSHF signifies, as in sect. 8, that both local- 
field effects and the electron-hole attraction have been included in the screening 
properties entering eq. (9.13). The Hartree-Fock band structure from ref. [49] is 
also shown for comparison. In fig. 7 the quasi-particle band structure TDSHF is 




Fig. 7. - Quasi-particle band structure of diamond calculated within the GW 
approximation with the screening taken within the TDSHF approximation (full line) and 
within the RPA approximation (dotted line) (from ref. [19]). 

compared with the RPA simplified form obtained by neglecting the electron-hole 
attraction in the screening properties. Several interesting features can be 
inferred from these calculations. 

i) Self-energy corrections to the Hartree-Fock eigenvalues are positive 
for holes and negative for electrons, thereby reducing the amount of energy 
required to produce an electron-hole pair. In particular, a drastic reduction of 
the diamond band gap at r occurs from the Hartree-Fock value (15.0 eV) to the 
final correlated value (7.4 eV) within the TDSHF approximation, in very good 
agreement with the experimental value. Notice also that the magnitude of self- 
energy corrections for the valence bands is larger than for the conduction bands, 
thereby matching a general «rule of thumb » that to larger effective masses 
correspond larger self-energy corrections. 

ii) The magnitude of the self-energy corrections within the RPA 
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approximation is smaller than the corresponding value within the TDSHF 
approximation; in particular, the RPA band gap is increased by about 1 eV 
compared with the full calculation. This is consistent with the ab initio values 
obtained for the long- wavelength static dieleqtric constant e within the RPA 
(4.25) and TDSHF (6.10) approximations, respectively, thereby fulfilling the 
other «rule of thumb» that self-energy corrections increase with e . Thus, 
although it formally violates conservation criteria, the TDSHF approximation 
gives in practice better results for the quasi-particle band structure and for the 
optical properties than the RPA approximation (cf. also sect. 8). 

iii) At any given k, the magnitude of self-energy corrections increases 
away from the gap region, resulting in a narrowing of the valence bandwidth. In 
particular, the Hartree-Fock bandwidth (29.0 eV) is reduced to 25.2 eV by the 
TDSHF calculation, a value which is rather close to the experimental value of 
24.2 eV. To assess the importance of the energy dependence of the self-energy 
operator across the valence bands, an additional calculation has been performed 
where the energy dependence of the matrix elements (9.13) has been (almost) 
neglected by taking their values at E = and cutting off the energy integration 
of eq. (9.11) past the plasma resonance. The result is that, while the energy gap 
is found practically unchanged with respect to the full calculation, the valence 
band width now comes out to be essentially unrenormalized with respect to the 
Hartree-Fock value. This suggests that an energy-independent approximation 
works well at energies about the gap, while the energy-dependence of the self- 
energy becomes progressively important as soon as we leave the gap region 
moving down through the valence bands. This finding corresponds to an 
increasing weight of the plasma resonance in the energy integration of eq. (9.11) 
with respect to the low-lying electron-hole continuum, as one proceeds deep in 
the valence bands. To make a further check on this point, the contribution of the 
plasma resonance to the full calculation has been neglected by terminating the 
energy integration before the onset of the plasma resonance. No variation has 
been consistently found on the gap. Similar results have been found for 
silicon [42], confirming that quasi-particle properties in the vicinity of the gap 
are essentially entirely determined by coupling to electron-hole excitations. 

iv) The self-energy corrections of figs. 6 and 7 show a weak, although 
noticeable, fe-dependence, being in general more pronounced away from the 
centre of the Brillouin zone. This is consistent with the short-range property of 
the mass operator discussed after eqs. (9.10) and indicates that self-energy 
corrections are essentially dominated by large albeit rigid shifts of the valence 
and conduction bands. 

A comment on the self-consistency of the single-particle Green's function is 
in order at this point. For the purpose of evaluating the matrix elements of the 
mass operator, the single-particle Green's function entering eq. (9.11) has been 
represented in the form (8.48) with the bands taken from a previous Xa. 
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calculation which is closer to the final result than the Hartree-Fock calculation. 
Once the new quasi-particle band structure has been obtained by solving eq. 
(9.9), however, the calculation should in principle be repeated using again for Gi 
an approximate expression of the type (8.48), until self-consistency is attained. 
While in the original calculation for diamond this kind of self-consistency has not 
been achieved, later work [44] has verified that iteration causes only a small 
further change for the spectrum which is whithin the estimated numerical error 
(0.2 eV). 

We finally discuss the decay of (valence) quasi-holes in diamond which can be 
calculated whenever the approximation for the mass operator is energy 
dependent. Within the GW approximation the mechanism responsible for this 
decay is the Coulomb interaction among the electrons which leads to radiation- 
less transitions to final states with more than one hole in the valence bands 
(intra-band multiple Auger processes). Clearly, these processes must be 
consistent with energy conservation and the decay rate must vanish below a 
threshold, rising sharply beyond it. For diamond the decay rate can be estimated 
to become noticeable about 10 eV below the top of the valence bands. In fig. 8 the 
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Fig. 8. - Imaginary part of quasi-hole levels in diamond along two different symmetry 
directions within the TDSHF approximation (from ref. [19]). o along A, + along A. 



imaginary part of the eigenvalues E n {k) of eq. (9.9) is shown versus the real part 
along two symmetry directions. (Recall that the corresponding lifetime of the 
approximate (N- l)-particles states is given by (2Im{E n {k)})~ 1 .) The characte- 
ristic value of about 1 eV for this quantity sets an intrinsic limit of consistency on 
the validity of the quasi-particle approximation which has extensively been used 
in our calculations. 

Some of above results can be understood, at least qualitatively, in an 
analytical way. To this end, we represent the (time-ordered) polarizability 
matrix xt entering eq. (9.12) in the form (5.9) (where A and B are now identified 
with the density operator), and we perform the energy integration in eq. (9.11) 
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by suitably closing the contour in the complex energy plane. We obtain for the 
diagonal matrix elements of the remainder M' 



d(s F - EfcO 6(s,fc' - £ F ) 



(9.14) drdr'uf(r)M'(r, r'; z k )u k {r') = 



E 2 



drdr' wjf (r) u k .{r) v(r - r')p s (r') 



where a> s are the excitation energies of sect. 5. For s k at the top of the valence 
bands the first term within brackets at the right-hand side of eq. (9.14) 
dominates over the second one, yielding a positive self-energy shift, while for s k 
at the bottom of the conduction bands is the second term that dominates over the 
first one, yielding a negative self-energy shift. Combination of the two effects 
results then in a net narrowing of the band gap compared with the Hartree-Fock 
value. However, if one neglected the energy dependence of the screened 
interaction altogether and took W'(r, r'; E) outside the energy integral in eq. 
(9.11), one would obtain a self-energy correction 

(9.15) Jdrdr'Mftr)M'(r, r <. Sk ) Uk ( r ')^ 

2 2 



-2 20(*F-e*O 



j drdr' uf(r) u k (r) v(r - r') Ps (r') 



which is manifestly positive for both electrons and holes [51]. The fact that some 
sort of energy dependence of the screened interaction should be retained is borne 
out, in particular, in the prescription for obtaining the COHSEX 
approximation [5], which is sometimes used for comparison with the full 
dynamical calculations [41, 44]. 

To get additional insight, we may further simplify eq. (9.14) by neglecting 
the local-field effects in the screening (although they are important for getting 
good quantitative results), and by approximating the low-lying electron-hole 
continuum and the plasma resonance (in which the energy range covered by the 
dielectric function is usually partitioned) by two modes with dispersion wi(g) and 
oj 2 (q), respectively. We thus take 

(9.16) e T - 1 (g;^) = l+ , , °\ . A9 + 



Mg)-iSf ^-Wfl)-iof' 

where the sign of the infinitesimal imaginary frequency corresponds to a time- 
ordered quantity, and the coefficients Q : and Q 2 are related to the plasma 
frequency by the sum rule (E.5) 



(9.17) 



Q\ +Q 2 2 = <4. 
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The self-energy shift now becomes 
(9.18) / drdr' uf(r)M'{r, r'; Sk )u k {r') = 

f 2 
J dr ut(r) exp [iq • r] u k {r) ■ 

e k - s k , + - iS s k - s k - - co s (q) + i& 

which shows at glance that the valence bandwidth As val is shrinked by the plasma 
resonance for Ae val < a> 2 (.q). A corresponding analysis can be done for the 
imaginary part of the self-energy shift. 



2 D 2 
% 2<» s (q) 



10. - Screening of static impurities. 

The mass operator considered in the previous section was assumed to 
possess the full symmetry of the crystal. In real materials this symmetry is often 
destroyed by the presence of impurities which may sit either at crystalline sites 
(substitutional impurities) or at interstitial sites (interstitial impurities). In both 
cases, the electronic density rearranges about the impurities which, therefore, 
test the screening properties of the crystal both locally and at large distances. 

In the simplest case of a single weak impurity the theory of linear response to 
an external (static) scalar potential can be applied. In fact, the change in the self- 
energy operator to lowest order in the (bare) impurity potential v im9 (r) is given 
by 

c ( SG~ 1 (1 2)\ 

(10.1) 2,0., 2) -2X1, 2) = - J d3l S(l, 2)*(1, 3) + gff (3 ) W 3 > = 

= - #(1, 2)v imp (l) + / d3r(l, 2; 3)^(3) , 

where S is the self-energy operator of the perfect crystal. Equation (10.1) has 
been obtained by combining Dyson's equation in the form (B.4) with the 
definition (C.6) for the reducible scalar vertex function. Notice that, upon 
insertion of eq. (10.1) into eq. (3.14), the first term at the right-hand side of eq. 
(10.1) cancels with the potential term. One may thus regard the second term at 
the right-hand side of eq. (10.1) as the effective change in the self-energy 
operator over and above the perfect crystal expression [18]. 

The screened impurity potential can be introduced at this point by recalling 
the definitions (3.17) for the irreducible scalar vertex function and (3.19) for the 
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inverse (longitudinal) dielectric matrix. One obtains 

! c $C~ 1 (1 2) 

(10.2) j Mm, 2; 3) v imp (3) = - J d3 ^ V mp (3) = 

= - / d34 5G p ( 1 3 ; 2) ^ ) ^ imp (4) = /d34f(l, 2; 3) e ~\S, 4) Vimp (4) = 

= Jd3Al, 2; 3)^(3), 
where use has been made of the «chain rule» (A.8) and 

(10.3) £ imp (l) = Jd2 e~\l, 2)v imp (2) = j dr 2 e-\r lt r 2 ; w = 0)v imp (r 2 ) 

is the (statically screened) impurity potential. Since at zero frequency the time- 
ordered inverse dielectric matrix in eq. (10.3) coincides with its retarded 
counterpart (cf. sect. 5), eq. (10.3) could have alternatively been obtained within 
linear response theory as the change in the Hartree potential resulting from the 
induced charge density about the (static) impurity [52]. Retaining only the 
effects of the impurity potential (10.3) corresponds to approximating 
f(l, 2; 3) = £(1, 2)£(1, 3) in eq. (10.2), i.e. to neglecting vertex corrections. It 
has been shown [18] that for a Coulombic impurity, whereby the Hartree 
potential (10.3) is 0{r~ l ) at large distance r from the impurity, vertex corrections 
produce terms of shorter range [0(r" 2 )]. Within the GW approximation used 
previously, this feature results from the converge properties of the single- 
particle Green's function far from the impurity centre. 

We do not enter here into the question whether vertex corrections within 
linear screening or nonlinear screening effects (due to higher terms in the 
expansion (10.1)) might be relevant to the calculation of specific impurity 
levels [53]. Rather we focus on the impurity potential (10.3) (or, equivalently, on 
the corresponding induced charge density) and take advantage of the knowledge 
of e~\r, r'; co) developed in sects. 8 and 9 to determine the size of the corrections 
introduced by local-field and excitonic effects. Details of the calculations are 
given in refs. [33] and [54] where embedding a pointlike impurity with positive 
unit charge both at substitutional and interstitial positions in crystalline diamond 
and silicon was considered. Different levels of approximation are again identified 
by the acronyms TDSHF, RPA, and RPA where the latter now means 
neglecting the electron-hole attraction and taking e -1 diagonal in the reciprocal 
lattice vectors. Figures 9 and 10 show the induced-charge density along two 
different directions for a substitutional impurity in diamond and silicon, 
respectively. These results demonstrate that in diamond both local field and 
excitonic effects introduce corrections of comparable magnitude. In addition, 
since these corrections are of the same magnitude of the isotropic RPA 
screening, both effects are crucial for a microscopic description of the impurity 
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screening. In silicon, on the other hand, there is a smaller difference between the 
RPA and the TDSHF approximations. Similar trends show up also at large 
distances from the impurity centre as it results, e.g., from the value of the total 
induced-charge density Q ind = (1 - e )/e . For silicon one gets the values e = 10.9 
in the RPA, e = 10.3 in the RPA, and e = 11.6 in the TDSHF, that give values 
of Quia within a 1% range. For diamond, there is about a 10% difference among 
the corresponding values of Q ind . 

The fact that the electronic charge in silicon is more easily polarizable than in 
diamond manifests also in the strong anysotropy of the screening charge 
distribution which is more pronounced in diamond than in silicon. By analysing 
the contour plots in selected crystal planes [54, 33], one finds that charge flux 
toward the impurity in diamond is preferentially directed along the bonds, while 
in silicon this asymmetry is less pronounced. 

Finally, a parallel study conducted for an interstitial impurity shows that in 
both materials the screening is locally more efficient for the substitutional than 
for the interstitial case [54, 33]. 

Although the screened impurity potential (10.3) considered in this section 
was strictly static, other self-energy terms beyond the Hartree potential entail 
in general dynamical screening effects through the frequency dependence of the 
dielectric matrix. On the other hand, in the theory of excitons which is 



t 1 I 1 r 




distance from point charge (Lattice constant) 
Fig. 9. - Induced charge density, p^d, about a positive substitutional impurity in 
diamond, along two different directions (from ref. [33]). For the m eaning of various 

approximations see text, a) (1.1. Indirection, b) (l.O.O)-direction; •■• RPA, RPA, 

TDSHF. 
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distance from point charge (Lattice constant) 



Fig. 10. - Same as in fig. 9 for silicon. 

considered in the next section (whereby bound states originate from the 
interaction of an excited electron with the hole left behind), the static linear- 
response screening of eq. (10.3) is replaced by an «effective» screening function 
that depends self-consistently on the excitonic binding energy. 

11. - Excitonic states. 

We have seen in sect. 5 from the Lehmann representation of the density- 
density correlation function that the poles of this function occur at the exact 
excitation energies of the subset of states which are coupled to the ground state 
through the density operator. For a crystalline semiconductor above the 
fundamental single-particle energy gap £ g these poles merge into a continuum 
and one may apply the methods of sect. 8 to calculate the strength of the 
absorption coefficient at any given frequency. It is below this edge, however, 
that the most striking two-particle effects over and above single-particle 
excitations take place, with the appearence of sharp peaks in the absorption 
coefficient that correspond to bound electron-hole pairs [55]. This new kind of 
states is associated with truly isolated poles of the density-density correlation 
function and is produced by the same many-body mechanism which is responsible 
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for the shift of the single-particle peaks in the continuum toward the low-energy 
side. In this section we describe an algorithm to locate these poles, which 
reduces the Bethe-Salpeter integral equation for the two-particle correlation 
function to an effective eigenvalue problem [12, 10, 8]. The need for this kind of 
algorithm arises also to circumvent the uncertainties in the single-particle band 
structure that might in practice be even larger than the spacing of the bound 
excitonic levels we are looking for. We shall, in fact, regard the underlying 
single-particle band structure as given (at worst, by fitting the energy gap to the 
experimental value) and concentrate on two-particle effects by calculating the 
excitonic binding energies. 

The derivation of the effective eigenvalue problem starts from the analysis of 
the linear combinations of the four time variables in the two-particle Green's 
function G 2 (l, 2; 1', 2') to identify its particle-hole portion for the relevant types 
of excited states. One may take, in particular, the symmetric combinations 

(11.1) T^^-ty, r 2 = t 2 ~t 2 ,, T=t 1 -f 

with t 1 = (1/2) + t v ) and f = (1/2) (£ 2 + k), as independent time variables, and 
conveniently group the 24 terms of the two-particle Green's function into 6 
classes with 4 terms each, as discussed in appendix G. The 8 terms that contain 
the phase factors exp [i(E s - E ) t 1 ] and exp [i(E s - E ) t 2 ] are then called the 
particle-hole Green's function. Upon setting t 2 . = £ 2 + + ) and Fourier 
transforming the variable t 2 , we obtain 

(11.2) L ph (x h x 2 ; x r , x 2 .\t u t r ; tu) = 




dt 2 exp [- ict>y L ph (x! t h x 2 £ 2 ; x v t r , x 2 - 1$) = 




exp 



.(Es-EoM 
1 2h — 



+ 



+ exp[- ico(f + |tti|/2)] - 
n 

y xf(*2, x 2 >; -S)j^%xi, x r ; tQ 
, h(o + (E s -E )-ir) 



exp 



.( E.-E )\z! 
1 2tT~ 



where rj^> + and the transform of the single-particle factor in eq. (4.1) has been 
dropped, since it vanishes provided &> is nonzero. 
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The Bethe-Salpeter equation (4.5) can be Fourier transformed accordingly. 
Assuming a particular (simple) pole at hw s = E s — E > to be isolated from the 
other singularities in the complex to plane, we may pick up the corresponding 
residue by integrating over an infinitesimal contour enclosing the pole. We 
obtain 

(11.3) x? h (x 1; x r ; n) exp [- ico s t l ] = 

= j d33'44'G 1 (l, 3)^(3', l')5(3, 4'; 3', A) x f{x„ x v ; n) exp[-to.fl, 

which can be regarded as a (homogeneous) bound-state equation for the right- 
hand particle-hole amplitudes (G.3a). 

Equation (11.3) holds for hoj s <s g , whereby the single-particle term at the 
right-hand side of eq. (4.5) gives no contribution and the pole at ha> s lies on the 
real energy axis. It can, therefore, be applied in the present form to determine 
binding energies of valence excitons. For core excitons, on the other hand, for 
which the core hole decays by Auger effect and the associated poles have finite 
imaginary widths, the procedure to obtain an effective eigenvalue equation 
requires a suitable generalization [56]. 

Equation (11.3) holds quite generally regardless of the specific form of the 
effective two-particle interaction E which embodies the dynamics of the bound 
electron-hole pair. Similarly to the calculation of the optical properties described 
in sect. 8, we take E of the form (cf. eqs. (4.4), (3.11), (3.12), and (8.44)) 

(11.4) 3(3, 4'; 3', 4) « - ihi(S, 3') «X4 + , 4') v(3, 4) + 

+ ih8@, 4) #3', 4')W(3 + , 3'); 

we now keep, however, the full time dependence of the dynamically screened 
interaction W without resorting to the static approximation (8.45). In this way, 
we can address studying the effects of dynamical screening on the excitonic 
binding energy Ef over and above the effective-mass limit in which Ef is much 
smaller than the energy gap e g and E reduces to a Coulomb potential modified by 
the (static) macroscopic dielectric constant e of the insulating medium [10]. 

We further represent the single-particle Green's functions entering eq. 
(11.3) by the quasi-particle approximation (8.48) that we known from sect. 9 to 
work well whenever the relevant single-particle excitations occur in the 
neighbourhood of the energy gap. We can then: i) project both sides of eq. (11.3) 
onto u c (xi)u$(xy), where c and v stand for conduction- and valence-band 
quantum numbers, respectively, and make use of the assumed orthonormality of 
the single-particle wave functions; ii) take the limit as ti = 0" and drop the 
common factor exp[— iu) s ti~\; iii) multiply both sides by [s c — e v — fauj. The 
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result is 

(11.5) [e c - s v - hco s ] j dxi dx y u*(x x ) u v (x v ) x? h (xi, x v ; -S) = 

= - 1 dx 2 dx 2 ' ufixi) u v {x 2 ) v(r 2 - r 2 .) xf(x 2 ., x z ; -S) + 

+ J dx 2 dx 2 . u*{x 2 ) u v (x 2 ') | exp [io> s t/2] W(x 2 , x 2 -; t + ) ■ 
• [e(-z-) exp [i-cujii] + 0(- t) exp [i-z-fe - ft&> s )/7i]] x? h (x 2 , *2-; f) . 

For the states of interest to us we can insert into eq. (11.5) an approximate 
expression for the right-hand particle-hole amplitudes which is obtained from the 
definition (G.3a) by expanding the field operators in terms of the single-particle 
wave functions u k (x) introduced above and retaining the same approximation 
that has led to the energy denominators of eq. (8.48). We thus set 

(11.6) x f(x, x'; t) = - exp [to s |r|/2] £ u c ,(x) u$(x') A S (V, c') • 

v'c' 

■ [0(t) exp [- is c > t] + 6(- t) exp [- ie B - t]] , 

where the expansion coefficients A s (v, c) can be expressed as transition matrix 
elements in terms of creation and destruction operators 

(11.7) A s (v, c) = (N\ata c \N, s) . 

The unknown excitation energies hm s and the associated expansion coeffi- 
cients A s (v, c) of eq. (11.6) are then determined as the eigenvalues and the 
associated eigenvectors of the eigenvalue problem (11.5). In fact, inserting eq. 
(11.6) into eq. (11.5) and performing the integration over i yield 

(11.8) [z c - Sv -ho> s ~\A s (v, c) + ^A s (v', c')\ (c, v'\v\v, c') - 

v'c' I 

-i ] 4^exp[-to£]<c, v'\W(u)\c', v)- 



h((o s — (x)) — (e c . — s v ) + iS h{do s + co) — (s c — £„■) + iS J 
(£— »0 + ), where we have introduced the notation 
(11.9a) (c, v'\v\v, c') = j dxdx' u?(x)u$(x')v(r- r')u v (x)u c (x') , 
(11.96) (c, v'\W(a>)\c', v) = j dxdx'u*(x)u$(x') W(x, x'; <o)u c .(x)u v (x') , 
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for the Coulomb matrix elements and their (dynamically screened) exchange 
counterparts (cf. also eq. (8.56)). 

Equation (11.8) constitutes a self-consistent eigenvalue problem whereby 
the Hamiltonian matrix and its eigenvectors A s depend on the eigenvalue ha> s , 
akin to the eigenvalue problem (9.2) for single-particle energy levels 
supplemented by the condition (9.6). Moreover, the Hamiltonian matrix in eq. 
(11.8) is Hermitian consistently with the assumption that the eigenvalue hoj s lies 
on the real energy axis below the threshold of the single-particle continuum. 

Equation (11.8) should be read as follows. The first term describes an 
uncorrelated electron-hole pair for which the two-particle interaction is 
effectively switched off. Interaction effects are contained in the two terms within 
braces. The first term is independent of the screening properties of the medium 
and, being of short-range character, does not affect the excitonic binding energy 
when the excitonic radius is relatively large. This term is responsible for the 
longitudinal-transverse splitting of excitonic states and induces a mixing of 
different spin-orbit eigenstates, to be discussed below. The second term within 
braces in eq. (11.8) bounds the electron-hole pair since at large electron-hole 
distances it reduces to an attractive Coulomb potential screened by e - 
Dynamical screening effects occur at shorter distances. To show this, we recast 
this term in a more meaningful form by making use of the relation between 
the time-ordered (T) and retarded (R) density-density correlation functions 
(cf. sect. 5) 

(11.10a) Re{ Xl (r, r'; w )} = Re{* R (r, r'; w )} , 

(11.106) Im{^ T (r, r'; w)} = sgn(w) Im{x R (r, r'; &>)} 

(&> real), of the property 

(11.11) Im{ ZR (r, r'; - *>)} = - Im{ XR (r, r'; w )} , 

and of the Kramers-Kronig relation 

(11-12) Re{x R (r, r'; oj)} =— da/ ; , 

~ J co — 0) 

where P means that the principal part is to be taken in the integration of the 
singularity. We obtain 

(11.13) i J" ~exp[-iojS](c,v'\W(oj)\c',v)- 



■ + - 



h(a> s — w) — (s c . — £„) + iS h(oj s + to) — (s c — £ t ,0 + iS J 
dxdx' u*(x)u c (x) j dr"e^l v Ar, r"; w s )v(r" - r')u v (x')u*(x') , 
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where we have introduced the effective inverse dielectric matrix 



(11.14) 



e£,„.(r, r'; co s ) = S(r - r') + J dr" v(r - r") / ^ 




The diagonal (c = c' and v = v') elements of expression (11.14) can be interpreted 
as describing the density response of the system with a hole in the state v probed 
by an electron excited into a conduction-band state c. Equation (11.14) has to be 
contrasted with the ordinary inverse dielectric matrix of linear response theory, 
namely, 



whereby an external (test) charge probes the density response of the system in 
the ground state. The effective and the ordinary dielectric matrices coincide, 
however, in the static limit when the differences hco s - (s c > — e„) and 
hco s - (e c - e„0 can be neglected in comparison with the characteristic energies in 
the loss matrix - Im {x R (r, r'; &>)}, which are at least of the order of the energy 
gap £ g . In other words, whenever E B « s g the effective potential at the right-hand 
side of eq. (11.13) reduces to a statically screened Coulomb potential with 
asymptotic form e 2 \r — r'| _1 1 for large values of |r — r'\. Deviations from static 
screening need otherwise to be considered. 

Numerical solution to the effective eigenvalue equation (11.8) could be 
tackled by converting it into a local-orbital representation along the approach 
described in subsect. 8.4. This procedure, which is appropriate to intermediate 
binding [57], is sketched in the appendix of ref. [56]. Here we restrict ourselves 
to considering the leading corrections to the effective-mass limit of eq. (11.8), 
whereby one follows the steps yielding the usual effective-mass equation for 
excitons [55], but retains the full effective inverse dielectric function and not just 
its static limit. 

For a crystalline semiconductor, a Bloch wave vector Q can be associated 
with the many-particle state \N, s) of eq. (11.7), wherein the single-particle 
labels v and c stand for the triad of a band index n, a Bloch wave vector k, and a 
spin index a. Translational symmetry also requires k c = k v + Q in eq. (11.7) (up to 
equivalence in reciprocal space) and identifies Q as a good quantum number for 
the eigenvectors of eq. (11.8). In the following, we shall adopt for simplicity the 
model of a semiconductor with simple (i.e., nondegenerate) valence and 
conductor bands. Consistently with the effective-mass limit, we assume 



(11.15) 
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A s (cr„, <j c ; k v , Q) to be sharply peaked about a particular value of k v (say k ) that 
coincides with an extremum (maximum) of the simple band e v (k v ). By the same 
token, k v + Q is confined to a small neighbourhood of an extremum (minimum) of 
s c (k c ) which we take at vanishing Q for the excitonic state to be excited by an 
electromagnetic field. 

The eigenvalue equation (11.8) can now be transformed into a differential 
equation by multiplying both sides by (2n)~ s exp [ik v ■ r] and integrating over the 
whole Ar„-space, and by expanding the single-particle band eigenvalues about the 
relevant extrema up to quadratic order. One obtains 



(11.16) 



2 \ m; 



to? 



« r) + 



+ 2 TTTTj exp [i(&„ - k ) • r] — -^A s (t v ; ^\ K) ^ 



(2^ 



dft^ 

(2*) 3 



'An 



\v\<j v k v , <j C 'k V ' + Q) + 



+ i f ^J^exip[-i(oS]{<j c k v , c v .k V '\W(co)\<j C 'k v ,, cr v k v ) 

J OK 



I 



1 



Ef+TT-Z (*V - *o) 2 + (*. - *o) 2 + - # 



- + 



\ 



2m? 



+ - 



2m* 
1 



£ * + 2ro* ( *" _ kof + 2m* ( *"' ~ * o) ~ ^ ~ **y 

= - Ef F s (d v , <i c ; r) , 

where m* and m? are the effective masses associated with the (isotropic) valence 
band maximum and conduction band minimum, respectively, 



(11.17) 



Ef = £ c (k ) - e v {k ) - hco s 



is the excitonic binding energy referred to the relevant single-particle energy 
gap, and 



f dA: 

(11.18) F s {a v , <7 C ; r) = — ^ exp [i(k v - k ) ■ r]A,(s v , <j c ; k v ) 



is the envelope function. Notice that the limit Q—> has been explicitly indicated 
for the Coulomb matrix elements in eq. (11.16), since they behave pathologically 
in the neighbourhood of Q = 0. To extract this behaviour, the Bloch functions 
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entering the Coulomb matrix elements can be represented in terms of Wannier 
functions (cf. eq. (8.51) where a single local orbital is now associated with each 
simple band), yielding 

(11.19a) (<? c k v + Q, « v ,kM°vK °c'K + Q)~$^X„ c ,^- l J w {Q), 

(11.196) JJQ) = 2 exp [iQ ■ t\ \ drdr' fc(r) $*(r)v(r - r') $*(r' - I) $ c {r' - D , 
i J 

where only dipole-like excitations on the same lattice site have been retained. 
Assuming further that dipole like excitations at i ^ do not overlap with the one 
in the central cell (/ = 0), we can expand the Coulomb integrals in eq. (11.196) 
into multipoles and retain only the leading dipole-dipole interaction terms: 

(11.20) J K {Q) - 1 drdr' &(r) #(r) v(r - r') tf(r') &(r') + 

+ 2 exphQ • I] , 

1*0 1*1 

where we have introduced the dipole moment of the excitation v— »e 
(11-21) d K = e\&rj>*{r)rUr). 

For a cubic material in the small Q limit the dipolar sum in eq. (11.20) acquires 
the form [58] 

1-1 1 oo\ V „ r •/-» I n [dy C ' d vc ~ 3(? • d*) (l • rf w< ;)] 

(11.22) X exp [iQ ■ I] — = 

_ 4ti 

3O 

which is a piecewice continuous at Q = and has a rapid angular variation about 
this point. 

Within the above approximations, the second term at the left-hand side of 
eq. (11.16) becomes 

(11.23) 2 f £rk exp [i(k v -k )-r]\^- s , <r, ; Av) JTQ ■ 

' Urn {<j c k v + Q, a* k v \v\ a v k v , a c . Av + Q ) — 



3(Q-tf*)«MJ 

e 2 



d* d 
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This term, being of short-range character, does not affect the excitonic binding 
energy when the excitonic radius is relatively large. This term, however, is 
responsible for the longitudinal- trans verse splitting of the excitonic states, since 
eq. (11.22) has different values when d vc is parallel or orthogonal to Q [59], and it 
induces a mixing of different spin-orbit eigenstates that leads to an inversion of 
the 2 : 1 intensity ratio [60]. 

The dynamically screened exchange matrix elements (11.96) that enter the 
last term at the left-hand side of eq. (11.16) may instead be approximated by 
neglecting both the variation of the periodic part of the Bloch functions about the 
extrema and the local field effects. One obtains 

(11.24) (, c k v , ^k v .\W(<o)\ 7 r k v )^i,^i,^~W{k v -k,; w ) , 
where 

(11.25) W(q; = v(q) e^Kq; <o) 



is the Fourier transform of the dynamically screened interaction obtained by 
neglecting off-diagonal screening. Alternative procedures may be followed at 
this point to evaluate the integrals over wavevector and frequency in the last 
term at the left-hand side of eq. (11.16), yielding, e.g., the Haken potential (cf. 
appendix H) or the asymptotic dynamically screened Coulomb potential. To 
obtain the latter, we neglect the dispersion of the single-particle energy levels in 
the energy denominators and make use of a simplified form of eqs. (11.13) and 
(11.14). The last term at the left-hand side of eq. (11.16) then becomes 



(11.26) 2 



dk v 



dk . 

exp [i(k r - k ) ■ r] I — -^A s (^,, <v; Av) - t^o • 
(2-) 3 J (2-) 3 



dha 
2- 



exp[- iajo] (7 c k„ v Av|W(w)| av^v, ? v k v ) 



1 



m+ £,(k l ,-k y- + £,(k v -k o) + H,-u 



-+ 



i 



2m* 



2m* 



dq 

(2n? 



exp[iq-r]v(q)e-\q; Ef) 



F s (<j„, <j c ; r) , 



where we have introduced the wave vector dependent effective inverse 
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dielectric function 



(11.27) 



e" 1 (g;^) = l-- f dfuo 



(-Dining; «>)} 



Dynamical effects thus modify the statically screened electron-hole attractive 
interaction, which is usually considered in the effective mass limit for 
excitons [55], by replacing it with the more attractive interaction within 
brackets at the right-hand side of eq. (11.26). Notice that the static limit can be 
recovered by setting Ef = in eq. (11.27). A plot of e (Ef) = [e'Kq = 0; E?)]~ 1 is 
shown in fig. 11 for the loss function - Im {e^Kq; w)} corresponding to silicon. 




0.1 0.2 0.3 0.4 0.5 

F B (eV) 



Fig. 11. - Effective dielectric function (silicon) for vanishing momentum transfer versus 
the excitonic binding energy Ef (adapted from ref. [56]). 



Inserting the results (11.23) and (11.26) into eq. (11.16), we obtain eventually 
the desired differential equation for the envelope function 



(U ' 28) -%[& + -3)+ r **'i» + 



+ S(r)Q limJ vc (Q)S^ c £ F,fa, ^; r)S, t 



Q^0 



- / exp [iq • r] v{q) e ~\q; Ef)F s (<r v , cr c ; r) = - EfF s (a v , cr c ; r) . 



(2^) : 



The spin structure of eq. (11.28) can be organized into a singlet solution, for 
which F s ( \ , | ; r) = F s ( [ , {; r), and a triplet solution, for which 
F s ( t , t ; r) = - F s ( i , i ; r). In accordance with the results of appendix D, only 
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the singlet solution is affected by the longitudinal-transverse splitting. (In 
particular, only the singlet transverse excitonic state can be excited optically.) 

Numerical solution to eq. (11.28) shows that the inclusion of dynamical 
screening effects within the effective mass limit leads to a negligible increase of 
the binding energy over the corresponding static value. However, in the case of 
core excitons a suitable modification of eq. (11.28) to include band structure 
effects (such as intervalley mixing and centrall-cell corrections) in an effective 
manner by modifying the screened Coulomb potential at short electron-hole 
separation, leads typically to a (20 -r- 30)% increase of the binding energy due to 
dynamical screening effects over and above static screening [56]. Moreover, 
deviations from static screening lead also to a narrowing of the Auger width 
when passing from the core hole to the core exciton [56], an effect which is a 
fingerprint of dynamical screening, since it cannot be obtained within the 
framework of a static theory. 
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Appendix A 



A functional derivative identity. 

Let F[T(x)] be a functional and let SF[T(x)] be its first variation when a 
small change ST(x) is made in the argument function T(x). The functional 
derivative 8F/8T(x) is then defined as that function of x for which 

(A. 1) 8F[T(x)] = F[T(x) + 8T(x)] - F[Hx)~] = f dy^-STiy) , 

J (y) 

to first order in 8T. This definition can be considered as a generalization of the 
ordinary (partial) derivatives of a function F({Xi}) whereby the variables {xj 
(i- 1, AO are specified at N points of a lattice: 

(A.2) F{{ Xi + SXi}) - F({Xi}) = £ ^-Sx, . 

i = i 

This remark implies that the usual rules for differentiation can also be derived 
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for functional derivatives. We consider some examples. 

(A - 4) U) + Fmz)l^^ 

S F[Y(x)] 



(A. 5) iii) 



1 


~SF[Y(x)] 


G[Y(x)f 


SY(y) 



SY(y) 



iv) Let F[Y(&; x, y] be the inverse of G[Y(£); x, 2/1 in the sense that the 
following identity is satisfied: 

(A.6) | dzF[Y(0; x, z] G[Y(0; z, y] = S(x - y) 

for any Y(S). Then 

v) Let ^[G[T(x); y]l Then 

^ ; SF(z) i ay SG(y)SY(x)' 

which is referred to as the «chain rule». 

We are now in a position to derive eq. (3.9) of the text. To this end, we recall 
that all the [/-dependence in the generalized single-particle Green's function 
(2.9) is explicit in the S factor. Varying the potential U by a small amount SU and 
keeping only first-order terms then yields (cf. eq. (A. 5)): 

(A* xm 9)- i <N\T[S§n»*X2W) rn ? AN\T[SS]\N) . 
(A.9) 2) - - - {m s]\N) ~ Gl(1 ' 2) (N\T[S]\N) " 

Since SS appears only within a time-ordered product wherein boson like 
operators, suchas Y\x, t + ) Y(x', t), can be treated as though they commute, we 
can evaluate SS as follows: 

(A.10) SS = - ^ | dtj dxdx' Y f (x, t + ) SU(x, x'; t) Y{x', t) . 

Entering eq. (A. 10) into eq. (A.9) and recalling the definition (2.10) of the 
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generalized two-particle Green's function result into 

(A.11) Sdd, 2) = | dtjdx dx' SU(x, x'; t) [- G 2 (l, x' t; 2, xt + ) + 

+ G { (l,2)G 1 (x't, xt + )]. 

Comparison with eq. (A.l) eventually identifies the expression within brackets 
in eq. (A.11) with the functional derivative of Gid, 2) with respect to U(x, x'; t), 
namely, 

(A. 12) 2) = - G 2 (l, x' t; 2, if) + G^l, 2) G,(x' t, xt+) . 

6U(X, X ', t) 

For a local potential of the type (2.5) we can write instead 

(A.13) ° Gl r!; Q f =-G 2 (l, 3; 2, 3^ + Gjd, 2)GA 3 + ). 

The generalized single-particle Green's function can thus be considered as a 
generating functional for higher-order Green's functions, while the external 
potential U acts as a source. 



Appendix B 



Dyson's equation. 

In the absence of both the interaction among the electrons (v) and the 
external potential ([/), the (zeroth-order) single-particle Green's function Gf 
satisfies the equation 



(B.l) 



Gfd, 2) : 



-ih^-~h{2) 

df 2 



Gfd, 2) = S(X, 2), 



whereby the corresponding inverse Green's function is given by 



(B.2) G{ 0) -\1, 2) : 



Oil 



<?(!, 2) = 



-ih^--h(2) 
dt 2 



8(1, 2) . 



The introduction of Gf ] 1 allows us to express eqs. (3.14) of the text as follows: 
(B.3a) { dS[G[°>-Kl, S)-U(l)$(l, 3) -2(1, 3)]Gi(3, 2) = i(l, 2), 
(B.36) Jd3Gi(l, 3)[G ( 1 0) - 1 (3, 2) - £7(2) £(3, 2)-S(S, 2)] = S(1, 2). 
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Equations (B.3a) and (B.36) identify the expressions within brackets as the left 
and right inverse of G u respectively. Demanding these two inverses to coincide 
leads thus to the identification 

(B.4) Gi\l, 2) = Gf 2) - U(1)S(1, 2) -2(1, 2) 

and 

(B.5) 2(1, 2) = 2(1, 2). 

The integral form of Dyson's equation (3.14) can readily be obtained by 
multiplying eq. (B.4) from the left by G?\3, 1) and from the right by G x (2, 4), and 
integrating over 1 and 2. The result acquires the familiar form 

(B.6) G.d, 2) = GfKl, 2) + 1 d34 Gf\l, 3) [C/(3) £(3, 4) + 2(3, 4)] G^A, 2) . 



Appendix C 



Reducible vs. irreducible parts of the correlation functions. 

a) Scalar and vector vertex functions. Generalized Green's functions of 
the form (2.9) and (2.10) can also be defined when the interaction Hamiltonian 
(2.4) is replaced by, e.g., the interaction Hamiltonian between a system of 
charged particles and the electromagnetic field: 

(C.l) H'(t) = q f dx 9 (r, t)p(x) - 1 f dxA(r, t) -j{x) + f djcA 2 (r, t)p(x) . 
1 c J 2mc J 

Here p(r, t) and A(r, t) are the scalar and vector potentials of the electro- 
magnetic field, respectively, q is the charge of a particle in the system, and p(x) 
and j(x) are the density and current operators 

(C.2) p(x) = f\x)nx), 



(C.3) j(x) = ^- [¥\x) \¥(x) - (VF + (x)) f(x)] . 

In what follows we shall omit the diamagnetic term of eq. (C.l) since its 
contribution vanishes in the limit A— >0 which will eventually be taken. 

Following the procedure outlined in appendix A one can readily show that 

(C5) _*_ [(T ,_ T , )i(1> 3;2> 3 , )W . 
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It is convenient at this point to introduce the (reducible) scalar and vector vertex 
functions according to the definitions 

(C.6) F(l,2;3) = SGl " 1(1 ' 2) 



(C.7) r(l, 2;3) = - 



8(9P(3)) ' 
SGTKl, 2) 



These functions can be related to the two-particle correlation function 
L(l, 3; 2, 3') by eqs. (C.4) and (C.5) and rule iv) of appendix A, to obtain 

(C.8) L(l, 3; 2, 3 + ) = |d45G 1 (l, 4)Gi(5, 2)7X4, 5; 3), 

(C.9) ^[(V 8 - V,.)L(1, 3; 2, 3')W = /d45G,(l, 4)^(5, 2)r(4, 5; 3). 

Equations (C.8) and (C.9) could alternatively be considered as definitions for F 
and T [16]. 

The integral equations satisfied by F and r can be obtained by combining 
eqs. (C.8) and (C.9) with the Bethe-Salpeter equation (4.5) for L. The result is 

(C. 10) F(l, 2; 3) = <y(l, 3) ${2, 3) + f d4567 ^ ( . 1 .' 2 _ ) . G 1 (4, 6) £(7, 5)7X6, 7; 3) 

J SGi(4, 5) 

and 

(C.ll) r(l, 2; 3) = yL[(V 8 -V 8 .)*(l, 3')*(2, 3)] s .. s + 

+ f d4567 ^ ( . 1 .' 2 - ) . G 1 (4, 6)Gx(7, 5)r(6, 7; 3). 
J oG 1 (4, 5) 

Equation (CIO) has to be compared with eq. (3.18) for the scalar irreducible 
vertex function F whose kernel is SM/SG } instead of SS/SG^. The difference 
between the two kernels is just the Coulomb term (cf. eq. (4.6)) 

(C.12) £ c (3, 5; 4, 6) = - ihS(S, 4)i(5, 6)v(3, 6) 

with respect to which F is defined to be irreducible. The irreducible counterpart 
of the vector vertex function can also be introduced through an equation of the 
form (C.9) (see below). 

For a system which is homogeneous in space and time the integral equations 
(CIO) and (C.ll) can be rewritten using the Fourier representations (D.13), 
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once the spin variables have been eliminated by summing over £ 3 
(C.13) r(p; q) = l + \ ^SXp, p>; q) G,{p' + q/2) G 1 (p' ~ q/2)r(p'; q) , 

(C. 14) r{p; q) = ^ + f ^-E s (p, p';q) + q/2) G x (p' - q/2) I\p'; q) . 

Notice how these equations involve only the effective interaction (D.80) for the 
singlet channel since the four-momentum q corresponds to a Coulomb line that 
carries no spin. 

b) The irreducible two-particle correlation function. By breaking the 
effective two-particle interaction (4.4) as the sum of the Coulomb term (C.12) 
and of the remainder 

(C.15) S R (3,5;4,6) = 8M(3 ' 4) 



SGi(6, 5) 



(cf. eq. (4.6)), one can introduce the irreducible part L of the two-particle 
correlation function as the solution to the integral equation 

(C.16) L(l, 2;V, 2') = 

= ^(1, 2')GA2, 1') + Jd3456G 1 (l, 3^(4, l')S B (3, 5; 4, 6)L(6, 2; 5, 2') 

with kernel S R . Once L is known, the full two-particle correlation function L can 
be obtained by solving the additional integral equation 

(C.17) L(l, 2; 1', 2') = 

= L(1, 2; 1', 2') + /d3456L(l, 4; 1', 3)S C (3, 5; 4, 6)L(6, 2; 5, 2'). 

The system of equations (C.16) and (C.17) is, in fact, equivalent to the original 
Bethe-Salpeter equation (4.5) for L. 

Several quantities can be expressed in terms of L. We give some examples. 

i) The irreducible polarizability (3.25) is given by 
(C18) zd, 2) = -ihL{\, 2;1 + , 2 + ), 

in analogy with eq. (4.13) for the full polarizability in terms of L. 

ii) The irreducible scalar vertex function (3.17) can be obtained from 

(C.19) L(l, 3; 2, 3 + ) = Jd45Gi(l, 4)Gi(5, 2)f(4, 5; 3) , 

in analogy with eq. (C.8) for the full-scalar vertex function in terms of L. 
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iii) The irreducible many-particle T-matrix can be defined as 

(C.20) L(l, 2; 1', 2') - dd, 2') d(2, 1') = 

= | d3456G 1 (l, 3)d(4, l')T(3, 5; 4, 6)^(6, 2')Gi(2, 5), 

in analogy with eqs. (4.9) and (4.12). Alternatively, f can be obtained as the 
solution to the integral equation 

(C.21) f(l, 2; V, 2') = 

= ~ R (1, 2; 1', 2') + / d3456~ R (l, 4; 1', 3)d(3, 6)d(5, 4)7(6, 2; 5, 2'), 

which is the analogue of eq. (4.11) for the full T-matrix T. T and f can further be 
related through the identity 

(C.22) T(l, 2; 1', 2') = f(l, 2; 1', 2')-ih \ d33T(l, 1'; 3)W3, 3')r(2, 2'; 3'). 

Equation (C.22) implies, in particular, that for the triplet channel T coincides 
with f, while for the singlet channel an additional contribution due to 
polarization is present (cf. also appendix D). 



c) Density-current and current-current correlation functions. Together 
with the (time-ordered) density-density correlation function 

(C.23) y.(l, 2)=-i(N\T[p'(l)p'(2)]\N) , 

n 

we can also consider the density-current and current-current correlation 
functions 

(C.24a) 2) = - 1 (N\T[p'(l)j'(2)]\N) , 

n 

(C.246) zd, 2) = - 1 (N\T[j'a)p'(2)]\N) , 

(C.24c) xd, 2) = - { (N\T[fa)f(2)]\N) , 

n 

where 

(C.25o) p'(\) = p{l)-{N\p{l)\N), 

(C.256) i'(D=j(l)-<iV|i(l)|iV>, 

are the density and current deviation operators, respectively. More general 
definitions of the type (3.22) may also be considered. 
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The correlation functions (C.23) and (C.24) can be expressed in terms of the 
(reducible) vertex functions as follows: 

(C.26a) *(!, 2) = - MAX, 2; 1 + , 2+) = - ih \ d45 G,{\, 4) Gi(5, 1)1X4, 5; 2) , 



(C.266) 2) = -ih 



-^(V 2 -V 2 .)L(1, 2; 1 + , 2') 

2im 



2'=2 + 



= -ih Jd45G 1 (l, 4)GA l)r(4, 5; 2), 



(C.26c) z (l, 2) = -ih 



2im 



(V!-V r )L(l, 2;1', 2 + ) 



r = i + 



: x(2, D = 



■ih\ d45G!(2, 4)^(5, 2)r(4, 5; 1), 



(C.26<Z) *(1, 2) = -ih 



l' = l + 
2' = 2 + 



= — ih 



^-(Vi-Vr) fd45G 1 (l, 4)^(5, l')r(4, 5; 2) 

2vm J 



r = i + 



= -ift Jd4567r (7, 6; 1)G X (6, 4)^(5, 7)IT4, 5; 3), 
where F is the noninteracting part of r (cf. eq. (C.ll)) 



(C.27) 



r (l, 2; 3) = ^[(V 8 - V».)*(l, 3')*(2, 3)] 3 . = 3 . 



Equation (C.26a) follows from eq. (C.8), while eqs. (C.266) to (C.26d) follow from 
eq. (C.9). Notice also that the four equations (C.26) have a similar structure and 
could accordingly be expressed as a single tensor equation of rank 4. 

The irreducible correlation functions associated with x, X> aR d X can De 
defined in analogy with eqs. (C.26), whereby the two-particle correlation 
function L is replaced by its irreducible part L. We have already seen, in fact, 
that the irreducible density-density correlation function can be expressed as (cf. 
eqs. (C.18) and (C.19)) 

(C.28) *(1, 2) = -ihL{l, 2; 1 + , 2 + ) = -ih Jd45G 1 (l, 4)^(5, l)r(4, 5; 2). 
We can similarly define 



(C.29a) 



~%X, 2) = -ih 



2im 



(V 2 -V 2 .)L(1, 2; 1 + , 2') 



2' = 2 + 
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(C.296) %l,2) = -ih 
(C.29c) %X,2) = -ih 

2' =2 + 

Furthermore, in analogy with eq. (C.28), we can also introduce the irreducible 
vector vertex function t (1, 2; 3) defined as the solution to the integral equation 
obtained from eq. (C.ll) by replacing the kernel SE/SG 1 with SM/SGi. In terms of 
r we can then express 

(C.30o) #1, 2) = %2, \) = -ih\ d45G,(l, 4)^(5, l)f(4, 5; 2) , 
(C.306) 3(1, 2) = - ih Jd4567r (7, 6; 1)^(6, 4)G X (5, 7)f(4, 5; 3) . 

Once the irreducible correlation functions (C.28) and (C.29) have been 
determined, their reducible counterparts can be obtained from eq. (3.24) and 
from 

(C.31o) xd, 2) = ^(1, 2) + Jd34z(l, 3)t*3, 4)$(4, 2), 

(C.316) zd, 2) = ^1, 2) + /d34Hl, 3)v(3, 4) #4, 2). 

Appendix D 

Elimination of the spin variables and transformation to the energy- 
momentum representation. 

The value of the spin variable <? does not change along a single fermion line 
whenever the Hamiltonian is spin independent. This implies, in particular, that 



(D.l) dd, 2) = 6 \ j2 G(x 1 , x 2 ), 

(D.2) M(l,2) = 6\, 2 M(x u x 2 ), 

(D.3) | d&r(l, 2; 3) = ^ & fX»i, x 2 , x») , 
(D.4) Jd&d&zd, 2) = x(x u x 2 ), 

(D.5) /d&e-Hl, 2) = €- 1 (x 1 , x 2 ), 



^-(V 1 -V 1 .)L(l,2;r i 2 + ) ii _ i+ , 



^-(Vi-V^^-^-V^Ld, 2; 1', 2') 

2im 2im r = i + 
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where we have introduced the notation x = (r, t). These identities allow us to 
readily remove the spin variables from the set of coupled equations discussed in 
sect. 3. Elimination of the spin variables from the Bethe-Salpeter equation, on 
the other hand, requires to distinguish between the two alternative paths that 
can be followed by joining the vertices of the effective interaction S. Specifically, 
we can set 

(D.6) 5(3, 5; 4, 6) = \ Si ^S\(x z , x 5 ; x 4 , x 6 )S kh + S hk S u (x 3 , x h ; x t , x 6 )& }(i , 

where s] and £ tt correspond to pairing the fermion lines in the way shown in 
fig. 12a) (the factor 1/2 in front of s\ has been introduced for later convenience). 
For example, the Hartree-Fock approximation for the self-energy yields 

(D.7) S H f(3, 5; 4, 6) = S hSi [- ihS(x s , xd v{x z , x 6 ) S(x 5 , x£)] + 

+ S isk [ihS(x 3 , x 6 )v(xi, xjSfa, a5 B )]& 4 { 5 , 

whereby the Hartree and the Fock terms correspond to a £} and to a 2 )t 
pairing, respectively (fig. 126). Alternatively, we could introduce the effective 
interactions for the singlet and the triplet channels 

(D.8a) S s ({xJ) = £ T t {{ Xi }) + S]({xi}) , 

(D.86) E t ({x i }) = S u {{x i }), 
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g. 12. - a) Alternative pairing of fermion lines in the effective interaction; b) lowest- 
order (Hartree and Fock) terms in the effective interaction; c) Bethe-Salpeter equation 
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and express the spin dependence of the (total) effective interaction S in terms of 
the Pauli matrices as follows: 

(D.9) E(S, 5; 4, 6) = ^S, 3 ; 4 d kk E s (x 3 , x 5 ; x 4 , x 6 ) + ^a h(i ■ <r kk S t (.x 3 , x 5 ; x 4 , x 6 ) . 

The T-matrix can similarly be decomposed. From the Bethe-Salpeter 
equation (4.11) one then finds that T u is obtained by iterating a T T only, while T\ 
derives from both E 1 and £ T t and couples also with T T T as shown graphically in 
fig. 12c). The coupled equations for T\ and !T T r can, however, be transformed by 
introducing the ^-matrices T a (a - s, t) for the singlet and triplet channels in 
analogy with eqs. (D.8), yielding two uncoupled equations for T s and T t : 

(D.10) T a (xi, x 2 ; x v , x 2 ) = Z a {x l , x 2 ; x v , x 2 ) + 

+ j dx 3 dx 4 dx 5 dx 6 £ a (x 1 , x 4 ; x v , x 3 )G x (x 3 , x^)G x (x b , x 4 )T a (x 6 , x 2 ; x 5 , x 2 ) . 

The transformation to the energy-momentum representation can be done 
after the elimination of the spin variables has been performed. To this end, we 
adopt the following convention for the Fourier representation: 

(D-ll) fix) = f-^L exp [ik ■ x] f(k) , 

where k = (k, co), dk = dk dco, and k ■ x = A; • r — cut. In particular, for any system 
which is translationally invariant in time and space we can write 

(D. 12a) M(x u x 2 ) = M{x x - x 2 ) , 

(D. 126) f(x u x 2 ; x 3 ) = f(x l - x 3 , x 2 - x 3 ) , 

(D.12c) ^(x 3 , x§, x 4 , Xg) = — (x 3 — Xg, X4 — Xg, X5 — x$) . 

This implies the following Fourier representations: 

(D. 13a) M(x u x 2 ) = f exp [iq ■ (x, - x 2 )] M(q) , 

(D. 136) r(x 1; x 2 ; x 3 ) = f -p- exp [i(p + q/2) ■ (x, - x 3 )] • 

• exp [- i(p - q/2) • (x 2 - x 3 )]f(p; q) , 

(D.13c) S{x 3 , x 5 ; x 4 , x 6 ) = f £ p £ q ' eXp[^(p+q/2)■X 3 ]eXp[^(p'-g/2)•X5]• 
-' (Z7T) (Z7t) (Z^f 

• exp [- i(p - q/2) ■ x 4 ] exp [- i{p' + q/2) ■ x 6 ]E(p, p'; q) , 
which are shown graphically in fig. 13. 
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Fig. 13. - Structure of M(q), f(p; q), and S(p, p'; q) from energy and momentum 
conservation. 



Elimination of the spin variables for the two-particle correlation function 
proceeds as follows. The functions L(xi<r u x 2 u 2 ; x v <j v , x 2I <s 2 ) corresponding to 
the six possible arrangements of the spin projections {a;} at the four external 
legs (cf. fig. 2a)) can be expressed as linear combinations of the singlet (s) and 
triplet (t) amplitudes 

(D.14) L S:t (x u x 2 ; x r , x 2 ) = L(x x f , x 2 f ; x v f , x 2 . f ) ± 

± L{x\ | , x 2 1 ; Xy f , x 2 ' I ) , 

with the + (-) sign referring to singlet (triplet), provided the ground state \N) 
possesses overall spin rotational symmetry (i.e., it is a spin singlet). In this case 
one may show that L s and L t satisfy two distinct Bethe-Salpeter equations with 
kernels S s and S t , respectively. 



Appendix E 



The Thomas-Reiche-Kuhn sum rule. 

In this appendix we show that a «$-derivable» approximation preserves the 
information about the (ground state) expectation value of the equal-time density- 
current commutator, thereby fulfilling the so-called Thomas-Reiche-Kuhn sum 
rule. In fact, a combination of eqs. (7.9) and (7.10) yields 

(E.l) A (N\T[p'{l) <F(3) y + (3')]|tf > + Vx • (N\T[j'(X) #(3) ^(3')]^) = 

= (*(!, 3') - 3)) (N\T[f(3) F + (3')]|iV) ; 
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applying then the operator (h/2im) (V 3 - V 3 to both sides of eq. (E.l) and 
setting 3' = 3 + leads eventually the equal-time density-current commutator in 
the form 

(E.2) (N\[p(x), j(x')]\N) = ^ (N\p(x')\N) T S(x, x') . 

The connection between the Thomas-Reiche-Kuhn sum rule and the commutator 
(E.2) can now be readily established by expressing (in the absence of magnetic 
fields): 

(E.3) Z(E s -E )\ Ps (k)\ 2 = 

s 

= ^ / dxdx' exp [- ik ■ r'] (N\[[ftx'), H\ p(x)]\N) exp [ik ■ r] = 
= || jdxdx' exp[ik-(r-r')]T-(N\[p(x),j(x')]\N), 

where V is the volume occupied by the system, E s -E are the excitation 
energies of the system, and 

(E.4) Ps (k) = V- 1,2 jdxex V [-ik-r](N\p(x)\N, s) 

are the corresponding matrix elements of the fluctuation density operator of 
wave vector A; between the ground state \N) and an excited state \N, s). 
Insertion of eq. (E.2) into eq. (E.3) gives finally 

(E.5) 2 (E s - E )\ Ps (k)\ 2 ="f^ y / dx(N\p(x)\N) , 

for any «$-derivable» approximation. Since the integral on the right-hand side of 
eq. (E.5) is taken to give the total number 2V of particles which are present in the 
volume V, eq. (E.5) reproduces the desired sum rule. In particular, taking 
k = (0, 0, k) with A; small, one can approximate 

(E.6) pM^-ikV-^q-'iNlllN, s) , 

where d z is the (^-component of the) dipole moment operator and q is the charge 
of a particle. Equation (E.5) thus acquires the more familiar form [61] 

(E.7) Z(E s -E )\(N\d,\N, s )f=^qK 



Well-known examples of conserving approximations for which the sum rule (E.7) 
can be explicitly verified are the RPA and the RPAE approximations [62]. 
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Appendix F 



The Clausius-Mossotti relation. 

The formalism presented in subsect. 8'4 provides a link between calculations 
based on the microscopic properties of the material and the macroscopic theory 
of the dielectric constant e M (w). One then expects to recapture from it the 
classical Clausius-Mossotti relation for the macroscopic dielectric constant 

1 +-~7ra(«) 

(F.D € M (o>)= | 

1 -— 7Wc(fc>) 

in the limit of completely localized electrons, such that the electronic excitations 
are represented by the polarization of point dipoles at the (Bravais) lattice sites 
(that we take throughout of cubic symmetry). In eq. (F.l) <x(w) is the site 
polarizability per unit volume that relates the polarization density (cf. eq. (8.8)) 
to the local electric field acting on the individual point dipoles [63]. The single 
function a(<u) thus provides the eigenfrequencies of the longitudinal and 
transverse excitations in the long-wavelength limit by considering the zeros of 
the numerator and of the denominator in eq. (F.l), respectively [64]. 

In the limit of extreme localization, we may retain only the form factors 
(8.53) with 1=0 and approximate 

(F.2) Mq + 0)~-iq-p e , 

where p c is the dipole matrix element 

(F.3) Pc = jdrtf(r)rh(r) 

between an occupied (v) and an unoccupied (^) local orbital. The net charge 
corresponding to this dipole excitation is thus zero. Notice that the integral in 
eq. (F.3) is limited in practice within the volume of the reference cell centred at 
r = 0. Combining eq. (F.2) with eqs. (8.26) and (8.59) yields for the macroscopic 
dielectric function 

(F.4) e M (a>) = 1 ZpI&Aq = 0; <*)pP , 

"0 cc' 

where /3 denotes any Cartesian component of the vector (F.3). The proper 
screening matrix S given by eq. (8.60) can be considerably simplified in the limit 
of extreme localization, by expanding the Coulomb matrix elements of eq. (8.61) 
into multipoles. To leading order, we get the dipole-dipole interaction 

(F.5) VUq) ~ \v%(q) = PtPc -q + 

+ r2 exp[-iq-m] ■ ; — tt , 
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where the prime on the summation means that the term m = is omitted. In 
other words, we have assumed that the exchange Coulomb matrix (8.56) (and 
possible modifications to it due to electronic correlation effects beyond the 
TDSHF approximation considered in subsect. 8"4) removes the self-interaction 
term from the proper Coulomb matrix (8.61) [65]. Proviso of this self- 
polarization correction, albeit irrelevant if the electronic wave functions are 
extended throughout the crystal, becomes indeed essential when the wave 
functions are localized [66]. 

The right-hand side of eq. (F.5) can be evaluated by the Ewald's method. 
Taking advantage of the results given in sect. 30 and appendix VI of ref. [28], we 
split the dipolar sum into a term that absorbs all that is nonregular at q = and a 
regular remainder. The nonregular term corresponds to the « macroscopic field » 
and cancels with the first term at the right-hand side of eq. (F.5), while the 
regular term corresponds to the so-called «inner field» and gives 

(F.6) VUQ = 0) - \vf,{q = 0) - ~^P? Pc ■ 

Insertion of eq. (F.6) into eq. (8.60) and of the result into eq. (F.4) requires us to 
sum the geometric series 



(F.7) |^t(? = 0>)#42^ 

"0 cc' "0 cc' 



(-ih)NUq = Q; «) 



2 (- ih)N°,Aq = 0; *)^-pf -pA~ ih)N° r Aq = 0; *>) + ... 

c"c" OLJq 



"0 cc- 



1 ~~ 77- 2pU- ih)N° c «Aq = 0; <o)p? + ... 



1 + ~3~X ^ 



where 

(F.8) x°(co)^^ r 2p^-ih)N° c Aq = 0; w )p^ = lim^x°(g + 0, q + 0; w ) 

is defined in terms of the polarizability for noninteracting particles (cf. eq. (8.58) 
with iV replaced by N°). Combining eq. (F.7) with eq. (F.4) leads eventually to 
the Clausius-Mossotti relation (F.l), provided we identify - e 2 /°(w) with the site 
polarizability per unit volume a(&>), computed without regard to the dipole 
interactions [64]. To this end, it is sufficient to consider the polarization density 
at the reference cell produced by the (self-consistent) local electric field 
(whereby the self-polarization correction has been subtracted [66]) set up by a 
purely longitudinal perturbing electric field. We thus write 



(F.9) P(r = 0; w ) = -J- / drr Pind (r; w ) = 

ao 
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= a~ S drr l dr '* 0(r - r ' ; w >^ (r '; w)= 

£i a 

= UrrUr'xXr, r'; o>) [^(r = 0; o>) - r' • £ loc (r = 0; «) + ...] = 

u q a 

= -4-1 drdr'rf (r, r'; w ) r' • # ]oc (r = 0; w ) + ... , 

where use has been made of the short-range property of z°(r, r'; &j) to confine 
both integration variables to the reference cell as well to neglect the spatial 
variation of the local electric field over the exent of the induced dipole, and of the 
cell neutrality condition [67] 

(F.10) {dr'z°(r, r';o>) = 

that follows directly from the orthogonality of the pair of localized orbitals 
entering eq. (F.3). The condition (F.10) can further be exploited to identify 

(F.ll) -±- \ArAr'rx\r,r';^)r'=x\^, 

thereby yielding the desired result 

(F. 12) P(r = 0; w ) = - e 2 %V) E loe (r = 0; «) . 

Appendix G 

The particle-hole Green's function. 

The two-particle Green's function 

(G.l) G 2 (l, 2; 1', 2') = ^- |J <iV|r[F(l) #(2) y + (2') F^l')]^) 

can be expressed as the sum of the following six functions: 
(G.2a) GKl, 2; V, 2') = |- lj <JV|T[VX1) ^d')] ZW2) F f (2')]|iV) • 

•^-|n-||t 2 | 

(G.26) G^l, 2; 1', 2') = |Y (N\T[Y(2) f\2')] T[¥(l) f\l')]\N) ■ 



1 



" T- 2 |Tl| ~2 |T21 
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(G.2c) Gl u (h 2; 1', 2') = 



■ ^ I <iV|2'tn2) Y\V)] T[T(1) T\2')]\N) ■ 



2 2 2 



fl T2 
- T + — + — 

2 2 



(G.2d) G 2 V (1, 2; 1', 2') = - 



- ^ (A^|T[F(1) T + (2')] 7W2) ^(l'XIlA 7 ) 



_ T2 _ J_ 

2 2 2 



2 2 



~1 ^2 
- + — + — 

2 2 



ft 



(G2e) G?(l, 2; 1', 2')= -f (A^riTd) F(2)] T[¥\2') Y\V)]\N) 



2 2 2 





^2 


1 




~2 














2 


2 


-7 + 


2 



ft 



(G.2/) g 2 vi (i, 2; r, 20 = -I (ArjTfip-^') nni) n2)]\N) 



2 2 2 



2 2 



Inserting a complete set of N-particle states (for functions I-IV), of (N + 2)- 
particle states (for function V), and of (AT - 2)-particle states (for function VI) 
between the two time-ordered products in eqs. (G.2), and introducing the right- 
and left-hand particle-hole amplitudes 

(G.3a) *f(x„ x,-; U - t 3 ) = (N\T[Y(i) f\j)]\ N, s) exp [i(E s - E ) (t t + t 3 )/2h] , 
(G.36) xf(x u Xj ; U - «,) = (N, s\T[Y(i) f*(j)]\N) exp [- HE, - E ) ft + ^/2ft] , 
the right- and left-hand particle-particle amplitudes 

(G.4o) x F(x t , Xf, U - tj) = (N\ T[f{i) Y(j)] \N + 2,s)- 

■ exp [i(E,(N + 2) - E ) ft + t,)/2h] , 

(G.46) x F(x t , x fi U - t,) = (N + 2, s\T[Y\i) f\j)]\N) ■ 

• exp [- i(E s (N + 2) - E ) ft + t,)l2h\ , 

and the right- and left-hand hole-hole amplitudes 

(G.5a) tfiXi, Xfi U - tj) = (N\T[Y\i) f\j)]\N - 2, s> • 

• exp [i(E s (N - 2) - E ) ft + t 3 )l2h\ , 

(G.56) ^ h (x„ Xj ; U - tj) = {N — 2, s\T[f{i) #(j)]|A0 • 

• exp [- i(E s (N - 2) - E ) ft + t 3 )l2h] , 
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(where the collective index s labels in each case the relevant set of states), we 
may rewrite the six functions (G.2) in the compact form 

(G.6) Gl(l, 2; v,2') = (-^X{a)6^-^\-^\T^- 

• 2 exp [i(E - E S (MS) tM ffcJgftTJ 

(a = I-VI). The entries needed in eq. (G.6) are listed in table I. 

Notice that the excitation energies E S (MJ - E enter eq. (G.6) only through 
the exponential factor that has been explicitly displayed. Suitable Fourier 
transformation with respect to the variable t, thus enables us to exhibit the poles 
of the spectrum E S (MJ - E . In particular, to locate the iV-particle excited states 
it is convenient to make use of the freedom in the choice of the time variables i 2 
and t 2 > in the Bethe-Salpeter equation (4.5) by setting £ 2 . = ^ + £ (£— »0 + ). With 
this choice and Fourier transforming the variable bi we may only retain the so- 
called particle-hole Green's function 

(g.7) Gta, 2; i', 2o=aa, 2; r, w+gihi, 2- v, 21, 

since the functions III- VI are nonsingular at the frequencies of interest. 



Appendix H 



The Haken potential. 

In this appendix we show how the often-quoted Haken potential, that 
contains corrective terms to the electron-hole attractive interaction - e 2 r -1 £q 1 
in the effective-mass approximation, can be derived from the last term at the 
left-hand side of eq. (11.16). To this end, we adopt the approximation (11.24) for 
the dynamically screened exchange matrix elements and take the (time-ordered) 
inverse dielectric function of eq. (11.25) within the so-called plasmon-pole 
approximation (cf. sect. 256) of ref. [5] and eq. (9.16) with a single mode) 



(H.l) efKq; *>) = ! + ' 1 ^ 



G„-l 



2 



where co = <o(q = 0). Equation (H.l) enables us to perform the frequency 
integration in eq. (11.16) in a close form. Changing in addition the integration 
variables k v and k v into p = k v — k and q-k v - k v> we obtain 

(H.2) 2 f i*£ exp [i(k v - k ) ■ r] f ~- B A s (^, *S, K) JVQ ■ 

■i exp [- i(oS\ (<r c k v , a v -k v . \ W(cu)|<j c .Av, <j v k v ) ■ 
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%o) + 2mt 



■ + 



i E *+^ - *o) 2 + ~^ (K - k f + hco - is 



+ ■ 



2raJ 



r dg r dp 

J \£lV) j \Lt7Z) 



dp 



(2*> 



x , ^o-l \ hcol 
€ ) 2w(g) 



-+ 



+ 



+- 



2tn? 2mjf 



Equation (H.2) may be further simplified by neglecting i) the p-dependence in 
the energy denominators owing to the localization of the coefficient A„ about 
p = 0, ii) the g-dependence of the dispersion relation &>(g) that we assume to be 
weak in comparison with the g-dependence of the single-particle energy levels 
about the extrema, and iii) the excitonic binding energy Ef in comparison with 
the characteristic screening energy ft&> - The result is the product of the envelope 
function F s (u V) <j c ; r) times the Haken potential (cf., e.g., ref. [59], sect. 76)) 



(H.3) v H (r)=-v(r) 



— + ^ 6 ° ^ (exp [- r/r e ] + exp [- r/r h ]) 



In eq. (H.3) r e and r h stand for the electron and hole polaron radii, respectively, 



(H.4) 



r h i 


1/2 




r h i 


1/2 


2m*&j 




2mja> 





and the quantity within brackets defines an r-dependent effective inverse 
dielectric function that interpolates smoothly between the values eo 1 for r-» °o 
and 1 for r = 0. It should be recalled, however, that the potential (H.3) is 
meaningful insofar as it describes the leading corrections for large r to the 
electron-hole interaction in the effective-mass approximation. 
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